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Finally,  differences  between  British  and  American  t hemal  data  for  nitro¬ 
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0.1$  in  estimated  nitrogen  contents  pointed  out  by  Raff  and  McHenry. 
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Introduction 


The  tables  originally  published  by  tbs  authorO)  in  1940  were  collected  from 
many  source s  and  had  quite  a  wide  distribution  in  this  country  owing  to  tbelr 
convenience.  A  similar  set  of  tables  was  published  independently  by  Hirsohfolder 
et  al.  in  America(2).  These  authors  included  corrections  for  departures  fran 
ideal  gas  conditions  at  the  pressures  obtaining  in  guns  but  such  corrections  have 
been  published  independently  by  Corner(3)  in  this  country  and  so  are  not  reproduced 
here.  The  properties  of  the  gaseous  pro  duo  ts  have  been  tabulated  at  100°  intervals 
and  smoothed  to  facilitate  interpolation.  Calculations  of  adiabatic  expansions 
are  considerably  simplified  by  the  use  of  a  table  of  entropies  and  so  this  has  been 
included.  A  table  of  free  energies  is  also  given  so  that  equilibrium  constants 
can  be  calculated  for  other  reactions  than  those  treated  here, 

information  on  sulphur,  potassium  and  chlorine  compounds  is  available  in  the 
open  literature  and  this  has  been  collected,  amplified  where  neoessary  and 
Systematised  in  a  GALCIT  Report  published  by  AllnanW.  in  independent  report 
by  Iuft(5)  also  includes  fluorine.  Consequently  no  information  on  these  compounds 
has  been  included  in  the  present  report. 

Some  of  the  data  for  the  product  gases  have  been  interpolated  from  that 
published  by  Wagman  et  al.(°) ,  some  has  been  obtained  by  modifying  older  published 
material  and  for  CO,  and  H,0  completely  new  calculations  from  spectroscopic  data 
have  been  carried  out.  full  details  are  given  in  the  appropriate  section. 

A  considerable  amount  of  information  on  heats  of  combustion  or  calorimetric 
values  of  propellant  ingredients  has  been  made  available  in  the  last  nine  years 
and  this  has  been  used  in  compiling  Table  1,  In  all  oases  care  has  been  taken 
to  give  as  accurate  a  value  of  the  calorimetric  value  as  possible  since  this  is 
very  nearly  equal  to  the  heat  liberated  in  a  propellant  explosion.  The  oal.val, 

(water  liquid)  is  a  very  convenient  quantity  to  tabulate  since  it  is  independent, 
to  within  1  oal./gm. ,  of  the  value  of  the  water  gas  equilibrium  constant.  It  is 
also  nearly  independent  of  the  temperature  at  which  the  measurement  is  carried  out, 
the  variation  being  of  the  order  of  -1  o&l»/gm.  per  15*0  increase  in  temperature. 

Corrections  to  oal.val,  determinations  have  been  discussed  in  a  previous 
paper(  7) ,  For  convenience  in  thermochemical  calculations  the  standard  end  state 
has  been  taken  to  be  CO* ,  CO,  H* ,  N2  and  water,  all  under  a  pressure  of  one 
atmosphere.  Published  values  of  cal.vals,  have  usually  been  corrected  for  the 
formation  of  byproducts  but  not  for  the  fact  that  the  final  cooled  products  are 
at  quite  a  high  pressure.  To  correct  to  one  atmosphere  pressure  we  mist  subtract 
a  quantity  proportional  to  the  density  of  loading  and  amounting  to  about  5  cal,/gpu 
for  a  density  of  loading  of  1  gm.  of  explosive  per  10  c.c,  bomb  volume,  this 
correction  varying  slightly  with  composition  of  the  explosive.  The  cal.vals. 
given  in  this  report  have  also  been  corrected  to  a  temperature  of  300*K,  The 
beats  of  formation  and  oombustion  in  Table  1  are  related  to  the  cal.vals.  by  the 
formulae 

,  4 

HJf.  =  67380  (0)  -  406 70  (C)  -  Oal.val.  (1) 

H.C.  =  134-720  (C)  ♦  33705  (H)  -  67380  (0)  *  Oal.val,  (2) 

these  figures  corresponding  to  heats  of  combustion  at  constant  volume  at  300*1  of 
94050,  674-10  and  67340  o&l ./mole  for  C,  Ba  and  CO,  The  corresponding  HJP.  of 
steam  at  constant  volume  at  300°K  is  57500  cal  ./mole.  These  are  conventional 
formulae  derived  as  follows  -  we  assume  (0)  >  (C)  so  that  ws  have  no  frpe  carbon 
but  2  (C)  +  £  (H)  >  (0)  so  that  wo  have  some  H*  and  CO,  If  no  COz  wore  formed 
we  should  have 

Cal.val.  =  26710  (C)  ♦  674-10  {(0)  -(C)]  -  H.F. 

while  if  no  H*0  were  forced  we  should  have 

Cal.val.  *  26710  (C)  ♦  67340  {(0)  -(C)]  -  H.F. 

Actually  we  get  rather  more  H*0  than  CO,  from  propellants.  If  we  assume  the  ratio 
of  the  two  amounts  to  be  4  :  3  than  our  formula  is  exaot,  but  the  error  is  always 
less  than  1  cal./gn. 
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It  should  he  noted  that  we  take  the  elementary  fora  of  carbon  to  he  graphite, 
the  other  elements  being,  of  course,  gaseous. 

The  oal.vals.  of  non-explosive  materials  such  as  oarbamite  are  their  effective 
values  when  small  quantities  are  added  to  a  propellant  and  in  many  oases  are 
negative  t  -  It  is  assumed  that  the  heat  of  miring  of  propellant  ingredients  is 
negligible, 

i 

1,  Heats  of  formation  of  propellant  constituents. 

(a)  Nitrocellulose .  3toe  original  tables  gave  figures  based  on  the  mean  values  of 
the  heats  of  combustion  given  by  Prettrev87  and  by  Schmidt  and  Becker (9),  m  eaoh 

case  samples  were  burnt  at  oonstant  volume  in  a  bomb  initially  filled  with  oxygen 
to  about  30  atmospheres  pressure,  Schmidt  and  Beoker  used  the  Lunge  method  for 
estimating  the  nitrogen  o  on  tent  of  the  N.C.  but  Prettre’s  method  is  not  stated, 
Schmidt  and  Becker  corrected  for  the  formation  of  nitric  acid  in  the  boob  and 
other  heat  corrections  are  negligible.  Apparently  they  did  not  correct  for  «h«lir 
in  the  nitrocellulose  which  might  amount  to  2$.  Such  an  amount  would  make  the 
observed  heat  of  combustion  too  low  by  1 2  oal./gn,  and  if  the  CO*  were  liberated 
in  the  nitrometer  it  would  add  .14$  to  the  measured  nitrogen  content,  Prettre's 
results  are  probably  subject  to  the  same  errors. 

The  results  of  both  sets  of  experiments  agree  quite  olosely  with  the  direct 
oal.val.  determination  of  de  Bruin  and  de  Pauw(l  0) ,  The  same  amount  of  ohalk 

would  here  produce  an  error  of  only  3  oal,/gm«  but  the  actual  challc  content  must 
have  been  less  than  this  since  the  N.C.  was  preoipitated.  These  authors  used 
the  Schulze -Thiemann  method  of  nitrogen  estimation,  which  aocording  to  Schmidt 
and  Becker  gives  results  some  .2$  lower  than  those  given  by  the  Lunge  method. 

The  oharges  were  fired  at  a  density  of  loading  "f  about  ,01  gn/c.c.  so  that 
corrections  for  byproducts  and  final  pressure  v.are  negligible;  the  bomb  was, 
however,  evacuated  to  2  ran.  air  pressure  before  firing  and  so  the  reaoticn  may 
have  been  inoomple  te . 

In  the  original  tables  the  mean  of  the  values  for  heats  of  oambustion  given 
by  Prettre  and  by  Schmidt  and  Becker  were  used  without  correction  and  in  oonsequenoe 
the  cal.val.  figures  deduced  were  probably  too  low.  More  recent  work  by  Milus(  1 1 ) 
gives  oal.vals.  of  N.C.  for  nitrogen  oontents  varying  from  12,62  to  13*45$.  The 
nitrogen  contents  mire  determined  by  the  du  Pont  nitrometer  method  and  so  are 
reasonably  accurate,  they  were  reproducible  to  ±  .02$  N.  No  mention  is  made  of 
chalk  or  ash  o  an  tent  but  since  the  measured  gas  volumes  agreed  very  olosely  with 
the  calculated  values  this  must  have  been  very  small.  The  density  of  loading 
was  0,1  gn/o.o.  and  subtracting  5  oal./gm.  to  bring  the  results  to  atmospherio 
pressure  we  get  by  least  squares 

Cal.val.  =  -  925  +  149.18  p  (3) 

where  p  is  the  $  by  weight  of  nitrogen  in  the  N.C.  This  fits  the  observed  results 
within  3  cal./gtn. 

Similar  work  has  been  carried  out  by  Taylor  and  Hall(l2)  for  four  nitrogen 
contents  ranging  from  12.1  to  13.24$.  In  three  cases  the  N,C.  was  preoipitated 
but  in  all  cases  the  ash  content  was  measured  and  allowed  for.  The  density  of 
loading  was  ,04  gn/c.o.  and  so  2  cal./gm.  must  be  subtracted  to  correct  the 
oal.vals.  to  abnospherio  pressure.  By  comparing  their  tables  9  and  10  we  see 
that  methane  formation  reduced  the  gas  volume  by  2.5  or  3  o.cs.  The  observed 
oal.vals.  were  therefore  too  high  by  3  or  4  oal./gn.  on  this  aocount.  A  further 
1  cal./gm.  must  be  subtracted  to  correct  to  a  temperature  of  27°C,  giving  final 
oal.vals.  of  1045  ,  974,  895  and  1054.  These  are  some  8  cal./gm.  lower  than  the 
figures  given  in  their  table  10  for  0°C.  Taylor  and  Hall  appear  to  have  increased 
their  oal.val.  figures  by  5  cal./gm.  to  correct  from  1 7°C  to  0°C  whereas  vre  should 
only  add  1  oal./gm.  and  so  it  seems  that  they  are  in  error.  Pitting  our  corrected 
figures  by  a  least  square  line  we  get  for  the  cal.val.  at  300°K  at  atmospheric 
pressure 

Cal.val.  =  -  799  +  140.06  p  (4) 

which  gives  the  final  figures  to  within  one  cal./ gn. 
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Measurements  of  the  heat  of  combustion  of  N.C.  have  also  been  made  by  Jessup 
and  ProsenO^ ).  These  authors  measured  the  heat  of  wetting  of  the  nitrocellulose 
which  decreased  from  4  cal./gm.  for  N.C.  of  11.5$  N  content  down  to  1  oal./gn,  for 
N.C.  of  13.4$  N  content.  Correction  to  atmospherio  pressure  and  correction  for 
nitric  acid  formed  reduced  the  observed  heat  of  combustion  by  3  to  4  oal./gm. 

The  impurity  present  in  the  N.C.  was  not  measured  directly  but  deduced  from  the 
ratio  of  the  mass  oi  CO2  formed  to  the  mass  calculated  for  pure  N.C. ,  the  deduoed 
amount  being  about  0,3$  for  N.C.  made  from  oottan  and  0.5$  for  N.C.  made  from 
wood  pulp.  The  nitrogen  determinations  were  made  by  the  suppliers  of  the  N.C. , 
the  Hercules  Powder  Company,  who  presumably  used  a  du  Pant  nitrometer  and,  as  is 
customary,  gave  their  nitrogen  oon tents  as  percentages  of  the  total  mass  of  N.C. 
plus  inpurities. 

For  pure  dry  N.C.  made  from  either  wood  or  ootton  the  heat  of  oombustion 
measurements  are  fitted  with  a  standard  deviation  of  less  than  1  cal./gm.  by 

H.C.  =  4179.50  -  140.7187  P  (5) 

and  so  the  cal.val.  is  given  by 

Cal.val.  =  -  806.76  +  139.725  P  (6) 

In  practice  the  N.C.  will  be  partly  wetted  but  the  heat  of  wetting  will  not 
usually  exceed  1  cal./gm.  of  finished  propellant.  Also  the  correction  to  27°C 

from  their  working  temperature  of  30°C  increases  the  heat  of  combustion  by  less 
than  3-  cal./gm.  In  using  equation  (6)  we  must  be  careful  to  use  the  true  nitrogen 
content  of  the  pure  nitrocellulose.  For  example  1  gm.  of  dry  N.C.  made  from  wood 
pulp  of  normal  standard  of  purity  and  of  nomina1  nitrogen  content  12.2$  will 
contain  about  .995  gms.  of  pure  N.C.  of  nitroge  1  content  12.26$.  If  the  impurity 
is  inert,  as  is  usually  the  case,  then  the  actual  cal.val.  will  be  902  cal./ gm. 
whereas  substituting  p  =  12.2  in  equation  (6)  would  give  898  cal./gm. 

Equations  (3),  (4)  and  (6)  give  the  following  values  for  pure  N.C. 


P  ($  N) 

11.5 

12 

12.5 

13 

13.5 

cal.val.  (3) 

791 

865 

940 

1014 

1089 

-  -  (4) 

812 

882 

952 

1022 

1092 

~  ~  (6) 

800 

870 

940 

1010 

1080 

One  would  expeot  Milas*  results  to  be  high  if  he  made  no  correotion  for  chalk 
and  other  inpurities  but  the  steady  difference  of  12  oal./gm.  between  Taylor  and 
Hall's  results  and  those  of  Jessup  and  Prosen  is  well  outside  the  experimental 
error.  If  we  assume  that  Jessup  and  Prosen  were  given  nitrocellulose 3  whose 
nitrogen  contents  were  assessed  on  the  pure  material  then  we  should  have  to  lower 
their  nitrogen  contents  and  increase  their  heats  of  oombustion,  the  nett  effeot 
being  a  reduotion  of  0.56  cal./gm.  in  the  heat  of  oombustion  or  cal.val.  for  a 
given  nitrogen  content.  Again,  ary  residual  solvent  would  increase  the  heat  of 
oombustion  and  reduce  the  cal.val.  Taylor  and  Hall  fired  their  charges  in  an 
evacuated  bomb  and  an  air  leakage  of  say  5  c.c.  per  gm.  of  ppt  would  not  only 
increase  the  measured  cal.val.  by  6  cal./gm.  but  would  also  increase  the  measured 
gas  volume  by  4  o.o./gm. ,  and  such  an  Increase  in  gas  volume,  which  would  have 
varied  considerably  from  round  to  round,  was  not  observed. 

The  discrepancy  can  be  explained  by  a  systematic  difference  of  about  .09 
between  the  percentage  nitrogen  contents  as  measured  in  the  U.S.  and  by  I.C.I. , 
the  latter  obtaining  the  lower  figure.  A  careful  comparison  of  the  two  methods 
by  Raff  and  McHenry (l 4)  showed  that  the  du  Pont  method,  as  used  in  the  U.S.  gave 
percentage  nitrogen  contents  consistently  higher  by  .10  than  those  obtained  by  the 
British  method.  The  du  Pont  method  avoids  a  number  of  systematic  errors  that 
arise  in  the  British  method  and  gives  the  theoretical  figure  for  inorganio 
substanoes.  We  shall  therefore  adopt  the  formula 

Cal.val.  =  -  810  +  140  p  (7) 

which  gives  the  same  values  as  (6)  over  the  range  of  interest,  and  make  the  proviso 
that  for  British  propellants  the  stated  percentage  nitrogen  content  should  be 
increased  by  0.10.  Equation  (7)  gives  cal.vals.  for  jure  N.C.  10  to  20  cal./gm. 


higher  than  those  in  the  original  tables.  The  figures  for  the  two  types  of 
nitrocellulose  oorananiy  used  in  British  propellants  with  nominal  percentage 
nitrogen  ocutents  of  12.2  and  13.15  are  912  and  1045  cal./gn.  respectively  and 
these  figures  oaribined  with  those  for  the  other  ingredients  give  a  nice  fit  with 
measured  oal.vals.  after  correction  for  byproducts  and  subtracting  the  pressure 
correction. 

To  avoid  confusion  we  mist  point  out  that  by  British  method  is  meant  the 
routine  method  used  in  British  propellant  factories,  which  is  intended  to  provide 
a  oheck  an  uniformity  of  production  rather  than  to  give  absolute  values  of 
nitrogen  content.  A  more  accurate  method  is  used  in  E.R.D.E.  for  special  purposes 
such  as  the  examination  of  nitrocellulose  reoovered  from  propellants. 

(b)  Nitroglycerine .  Sohmidt(l5)  and  RinkerihachO  6)  gave  1639  and  163O  cal./gm. 

for  the  heat  of  oanbustion  of  N.G.  at  ocnstant  volume  and  in  the  original  tables 
the  mean  of  these  two  was  taken  to  give  a  oal.val.  of  1784  cal./gm.  More  recent 
work  by  Burlot(17)  gave  a  heat  of  oanbustion  of  1623  oal./gn.  Unpublished  work 
in  the  Research  Department  by  Robertson,  using  mixtures  of  N.G.  and  gunoottcn  in 
a  propellant  calorimeter  gave  unoarrected  oal.vals.  for  N.G.  alone  of  1750  cal./gm. 
and  similar  work  by  Goodwin  and  Taylor  gave  mean  values  of  1 770  and  1 750.  .  Further 
unpublished  work  by  Sutton  and  Ambler  in  the  Research  Department  in  1932  gave  the 


following  results  for  the  oal.val.  of  N.G. 

N.G.  plus  stick  cordite  in  ppt.  bomb  1709 

N.G.  plus  ground  ........  1763 

N.G.  plus  Hs  in  heat  of  combustion  boob  1755 

From  oal.val.  of  60/40  ball istite  1738 

From  heat  of  combustion  of  60/40  ballistite  1770 


The  figure  adopted  by  the  calorimetric  section  for  calculating  the  oal.vals.  of 
propellants  containing  N.G.  was  1740  oal./gm.  the  corresponding  figures  for  N.C. 
(12.2^5  N)  and  guncotton  (13*1$  N)  being  923  and  1040  oal./gm.  respectively. 

The  most  recent  determinations  are  heat  of  combustion  measurements  by  Taylor 
and  HallO  2)#  Their  results  for  N.G.  were  not  so  regular  as  for  N.C. ,  probably 
on  account  of  oxidation  of  bomb  metal  by  the  excess  oxygen.  Their  best  estimate 
was  1615  oal./gn.  for  heat  of  combustion  which  gives  1761  cal,/ gm.  for  the  •oal.val. , 
corrected  to  atmospheric  pressure.  Their  corresponding  figure  obtained  by  firing 
propellant  mixtures  was  1750  cal./gn.  Since  the  combination  of  N.G.  and  N.C.  in 
a  gelatinised  propellant  is  probably  slightly  exothermic,  the  effective  value  of 
the  cal.val.  of  N.G.  will  be  a  little  less  than  1 761  cal./gm.  We  shall  therefore  ‘ 
assume  that  in  gelatinised  propellants  the  effeotive  cal.val.  of  N.G.  is  1750 
oal./gm. 

(c)  Diothylene  glycol  dinitrate.  The  heat  of  combustion  reported  by  RinkenbachO ®) 
was  2*798  oal./gm.  at  constant  volume ,  giving  a  cal.val.  of  1080  cal./gm.  Direct 
measurement  of  the  oal.vals.  of  nitrone  thane  and  a  60/40  nitrone  thane/D.G.N. 
mixture  gave  figures  of  1121  and  1079  oal./gm.  respectively.  The  heat  of  mixing 
was  -1  oal./gm.  giving  1051  cal./gm.  for  the  oal.val.  of  pure  D.G.N.  and  we  shall 
adopt  this  value. 

(d)  Dinltroglyool .  RinkenbachO £)  gives  1 764  cal./gm.  for  the  heat  of  oanbustion, 

which  equals  the  cal.val.  for  this  substance.  The  heat  of  combustion  given  by 
Tomioka  and  Takahashi(19)  is  1749  oal./gm.  The  oal.val.  of  a  65/35  mixture  of 
D.N.G./nitropropane  was  determined  in  the  A.R.D.  to  be  1021  cal./gm.  giving  a 
oal.val.  for  D.N.G.  of  1 730  oal./gm,  and  ws  take  1 740  cal./gm.  as  a  weighted  mean. 

(e)  Nitroguanidine .  guanidine  nitrate  and  oxamide.  Sehmidt05)  quotes  Matignan' s 
figures  of  2032  and  1715  cal./gm.  for  heats  o 7  combustionv  19)  of  the  first  two, 
leading  to  oal.vals.  of  737  and.  610  cal./gm.  respectively.  The  figure  for 
picrite  used  in  the  A.R.D.  calorimetric  section  is  715  cal./gm.  Other  heats  of 
oanbustion  quoted  by  Landolt  Bernstein  (e.g.  Ill  o)  for  picrite  are  2010  oal./gm. 
by  Tomioka  and  Takahashi  and  2018  cal./gm.  by  Badoohe. 

Some  very  precise  measurements  of  the  heat  of  oanbustion  of  Welland  picrite 
were  carried  out  by  Prosen  and  Gilmcmt(20)  who  obtained  2014  cal./gm.  Two  other 
samples  fired  by  Stegeman(2^  )  gave  heats  of  combustion  of  2002  and  2008  oal./gm. 


V 


respectively.  The  mass  of  CO2  formed  in  Prosen  and  Gilmcnt'a  experiments  was 
99.96#  of  that  calculated  for  pure  picrite  hut  the  actual  purity  of  the  sample 
was  not  known.  The  wartime  specification  for  hulk  production  allowed  about  1# 
of  impurities.  If  we  assume  that  the  only  impurity  in  Prosen  and  Gilmont'a 
sample  was  .C 4#  of  inert  matter  then  the  calculated  oal.val.  is  720  cal  ./gnu 
For  picrite  of  normal  purity  it  would  he  better  to  take  a  cal.val.  of  715  cal./gm. 
and  this  we  have  done. 

No  other  information  appears  to  he  available  for  guanidine  nitrate  and  so  we 
shall  assume  a  cal.val.  of  610  cal./gm.  The  only  available  figure  for  oxamide  is 
a  heat  of  combustion  of  2316  cal./gm.  given  by  Stohmann  and  HauemannC 22 )  giving  a 
cal.val.  of  -74-5  oal„/gn. 

(f)  Dlnltro toluene  and  Trinitrotoluene.  Schmidt05)  gives  4709  and  3599  cal./gm. 
for  the  beats  of  combustion.  Prosen  and  GilmontC2®)  were  supplied  with  D.N.T.  by 
the  Hercules  Powder  Company  which  produced  99.99#  of  the  calculated  mass  of  CO*  on 
combust  ion,  with  a  heat  of  combustion  of  4658  cal./gm.  StegemanC^I  > found  a  beat 
of  comhistion  of  3594  cal./ gm.  from  a  practice  firing  of  T.N.T.  of  unknown  purity. 
Heats  of  combustion  of  D.N.T.  and  T.N.T.  quoted  by  Tomioka  and  Takahashi  are  4691 
and  3592  cal«/gn.  while  BurlotO?)  gives  4660  and  3596  cal./gn.  respectively. 

Since  in  recent  years  only  T.N.T.  of  Grade  1  purity  has  been  used  old  figures  are 
probably  not  reliable.  Garner  and  Abernethyi22)  measured  the  heats  of  combustion 
of  all  the  isomerio  forms  of  mono-,  di-  and  trinitrotoluene,  their  figures  for 
2:4  D.N.T.  and  2:4:6  T.N.T.  being  4691  and  3623  cal./gm.  Their  figures  ranged 
from  4691  to  4730  and  from  3623  to  3677  cal./gm.  for  the  various  isomers. 

Prosen  and  Gilmont's  and  Garner  and  Abernathy' a  heats  of  combustion  of  D.N.T. 
give  cal. vale,  of  -150  and  -117  respectively.  The  figure  determined  by  Taylor 
and  Hall  from  the  oal.val.  of  propellant  mixtures  was  -135  cal./gm.  corrected  to 
atmospheric  pressure  and  we  shall  adopt  this  figure.  A  figure  of  -180  cal./gm. 
used  in  the  A.R.D.  calorimetric  section  seems  to  he  too  large.  Garner  and 
Abernathy's  figure  for  2:4:6  T.N.T.  gives  a  oal.val.  of  509  cal./gm.  Allowing 
for  impurities  in  grade  1  T.N.T.  it  is  best  to  take  a  figure  of  500  cal./gm. 

(g)  Diethyl.  Dibutyl  and  Dianyl  phthalate.  Ambler(23)  obtained  7389  and  7705 
cal./gn.  for  the  beats  of  combustion  of  dibutyl  and  dianyl  phthalate.  Schmidt's 
figures  were  a  little  higher  and  the  two  figures  were  interchanged  by  mistake. 

Pro  sen  and  GilmontC2^)  obtained  6419  and  7352  calVgm.  for  samples  of  diethyl  and 
dibutyl  phthalate  whose  carbon  contents  were  .9996  and  .9985  times  the  calculated 
values  for  pure  substances.  Taking  the  mean  of  the  two  figures  for  dibutyl- 
phthalate  the  corresponding  cal. veils,  of  the  three  substances  are  -1766,  -2070 
and  -2190.  Taylor  and  Hall  obtained  a  oal.val.  of  -2025  oal./gn.  for  dibutyl- 
phthalate  but  this  was  obtained  by  firing  a  sample  containing  only  2#  of  thia 
liquid.  These  large  (negative)  figures  are  not  likely  to  be  increased  greatly 
by  inpurities  in  the  materials  and  so  since  we  are  working  to  the  nearest  5  cal./ 
gm.  we  shall  take  -1765,  -2070  and  -2190  oal^/gm.  respectively  for  the  oal.vals. 
of  these  liquids. 

(h)  Stabilisers  and  Graphite.  Carbamite  is  the  British  name  for  die thyldiphenyl 
urea,  sometimes  called  ethyl  centralite  or  central! te  I.  its  heat  of  comtuetion 
is  given  by  Schmidt05)  as  8528  oal«/gn«,  110  oal./gm.  higher  than  Rirlot's 
figure  (17)  and  giving  a  negative  heat  of  formation.  Prosen  and  Gilmont  obtained 
8374  oal./gn. ,  leading  to  a  oal.val.  of  -2422  oal./gm.  This  is  in  fair  agreement 
with  Taylor  and  Hall's  cal.val.  figures  of  -2397  cal./gn.  derived  from  heat  of 
combustion  and  -2440  cal./gm.  obtained  from  the  oal.vals.  of  propellant  mixtures. 
The  last  is  the  moot  reliable  figure  and  we  shall  adopt  it. 

Tbs  oal.val.  of  methyl  central ite  was  determined  by  Taylor  and  Hall  as 
-2380  oal./gm.  The  value  deduced  from  their  heat  of  combustion  measurement  is 
-2303  cal./gn.  but  the  directly  determined  oal.val.  figure  is  the  one  that  we 
should  use.  Schmidt's  figure  for  the  heat  of  combustion  of  akardite  i.e. 
asymmetric  diphenylurea,  is  7700  oal./gn.  giving  a  cal.val.  of  -2140  cal./gn.  and 
a  negative  heat  of  formation. 

The  composition  of  mineral  jelly  is  not  known  accurately  but  the  C:H  ratio 
is  a  little  higher  than  that  calculated  for  a  heavy  aliphatic  hydrocarbon.  Some 
of  the  carbon-carbon  bands  are  therefore  double.  If  1#  of  the  carbon  bonds  are 
double  then  the  H  content  is  1#  less  than  for  a  saturated  aliphatic  hydrocarbon 
wwfl  the  negative  cal.val.  80  cal./gn.  smaller.  The  special  mineral  jelly  used 
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for  M.D.  cordite  had  a  bromine  number  of  70  to  110  i.e.  about  ,001  gn.  moles  of 
bromine  reacted  with  1  gn.  of  M.J.  This  implies  that  about  0.6£  of  tbs  C-C  bonds 
were  double.  For  propellants  containing  oarbamite  ordinary  mineral  jelly  of  U,S, 
origin  is  used,  for  which  the  bromine  number  is  given  by  Redwood(25)  as  about  10. 
We  may  therefore  take  present  day  M.J,  to  be  practically  pure  paraffin  containing 
15  to  20  carbon  atoms  per  molecule.  .  if  we  take  mineral  jelly  to  be  pure 
ootadecane,  Cte  Hse>  "wran  Pro  sen  and  RossiniC2®)  give  a  heat  of  combustion  at 
constant  volume  of  11256  cal./gn.  leading  to  a  cal.val.  of  -5505  oal,/gn.  For 
special  M.J.  used  in  M.D.  we  shall  take  a  cal.val.  of  -5260  cal./ gn.  and  a 
corresponding  reduction  in  fydrogen  content. 


For  diphenylamine  the  heat  of  combustion  is  given  by  Schmid tO  5)  as  9086 
cal./gn.  Figures  of  9081  ..9085  and  90^9  cal./gn.  are  given  by  Burlot\l7), 
Stohmann(27)  and  Lemoult(28)..  Kharasoh*s  calculated  figureC2?)  is  91 1 8 .  oal./gn. 
No  more  modern  figure  has  been  found  and  so  we  shall  take  the  mean,  giving  a 
cal.val.  of  -2660  cal./gn.  For  camphor  Kharasch  calculates  a  heat  of  combustion 
of  9285  and  quotes  an  unpublished  figure  of  927k  cal^/gn.  by  Swie to slawski . 

Taking  9276  cal./gm.  we  calculate  a  cal.val.  of  -26 70  cal./gn.  Schmidt(9)(l5) 


gives  heats  of  combustion  for  a  number  of  other  stabilisers  which  have  been  used 
in  German  propellants  but  which  are  of  little  inmediate  interest. 


The  cal.val.  of  pure  graphite  would  be  -5586  cal./gn.  but  the  graphite  used 
in  propellants  may  contain  up  to  1(^5  of  impurities.  We  shall  therefore  take  the 
figure  of  -5565  cal./ gn.  obtained  directly  by  Taylor  and  Hall. 


(i)  Solvents  and  Water.  The  cal.val.  of  water  is  zero  and  the  heat  of  wetting 
of  N.C.  in  propellants  is  negligibly  small.  The  heat  of  wetting  of  N.C.  by 
ace  tone  (3®)  less  than  1  cal./gm.  when  the  ace + one  weight  is  less  than  3%  of 

the  N.C.  weight.  Heats  of  wetting  by  ethyl  ethar  and  alcohol(31 )  are  2  cal./gn. 
for  and  3%  solvent  respectively.  All  these  figures  are  negligible,  for  the 
amounts  of  solvent  usually  to  be  found  in  propellants. 


.  The  heat  of  combustion  of  acetone  is  7282  7551  or  7M9  cal./gm.  as  calculated 

by  KharaschC2?),  measured  by  Berthelot  and  Delepine(32)  or  measured  by  finery  and 
Benediot(35).  Taking  the  mean  we  get  a  cal.val.  of  -1 950  cal./gn.  with  an 
uncertainly  of  50  cal./gm.  Fortunately  acetone  only  occurs  in  small  quantities 
in  propellants. 


Diethyl  ether  has  a  latent  heat  of  vaporisation  of  87  cal./gn.  at  20°C. 

Thomsen (3k;  and  Stohmann,  Rodatz  and  Herzberg(35)  measured  the  heat  of  combustion 
of  the  vapour  and  the  latter  also  measured  the  heat  of  combustion  of  the  liquid. 
Their  figures  corrected  to  the  liquid  state  at  constant  volume,  are  8811,  8818 
and  8789  cal./gn.  respectively.  Taking  the  mean  we  get  a  cal.val.  of  -2100  cal./ 
gn.  For  ethyl  alcohol  we  have  the  very  accurate  heat  of  combustion,  given  by 
Rossini (36^  of  7078  cal./gm.  corrected  to  constant  volume  at  25°C.  This  gives 
a  cal.val.  of  -1700  cal./gm.  in  round  figures. 

(  j)  High  explosives.  We  use  Schmidt*  s  heat  of  combustionC  29 )  of  PETN  which  agrees 
closely  with  Burlot*s  figure.  The  A.R.E.  Calorimetric  section  have  found  a 
cal.val.  of  1560  cal./gn.  for  RDX  in  propellant  mixtures  and  this  is  more  reliable 
than  figures  of  1575  cal./gm.  or  1570  cal./gn.  calculated  from  Schmidt*  s  or 
Bur  lot*  s  heat  of  combustion.  For  tetryl  both  Schmidt  and  Stegeman(21)  give  a 
heat  of  combustion  of  2925  cal./gn.  agreeing  closely  with  Bur  lot  and  leading  to 
a  cal.val.  of  950  cal./gn.  Stegeman  gives  heats  of  combustion  of  some  thirty 
other  explosives,  mostly  of  only  academio  interest.  A  comprehensive  collection 
of  information  on  explosives  is  given  by  Blatt(^-O). 

(k)  Liquid  propellants.  The  figures  for  methyl  and  ethyl  nitrate,  nitrome thane, 
nitropropane  and  8 5$ "hydrogen  peroxide  sure  taken  from  Wheeler  et  al.(57),  those 
for  nitric  acid  and  liquid  oxygen  from  Bichowsky  and  Rossini(58)  and  that  for 
decane  from  the  heat  of  combustion  given  by  Rossini(59).  The  figure  for  liqiid 
oxygen  is  the  heat  liberated  in  cooling  oxygen  gas  from  500°K  to  90°K  and  then 
liquifying  it.  Figures  for  hydros©11  peroxide  of  any  strength  and  for  calcium 
permanganate  axe  given  in  an  unpublished  report  by  Pike  and  Green. 

The  heats  of  formation  have  been  given  for  constant  pressure  as  well  as 
constant  volume  conditions,  for  convenience  in  rocket  calculations.  The 
difference  consists  in  adding 
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*RT  {(H)  +  (N)  +  (0)] 

=  298  {(H)  +  (N)  +  (0)|  cal./ga. 

Finally  the  heat  of  formation  of  diamond  is  given  for  use  in  correcting  heats  of 
formation  given  by  BichowBiy  and  Rossini(38). 

2,  Thermochemioal  data  for  the  products  of  explosion, 

(a)  Heats  of  formation.  The  heats  of  formation  of  CO*,  CO,  H*0  and  CH4  are  taken 
from  Wagnan  et  al.lt>),  that  of  OH  is  derived  from  Dwyer* s  data(41),  that  of  NO 
from  Biahowsky  and  Rossini(39)  and  that  of  NHj  from  Gillespie  and  Beattie(42). 
Heats  of  dissociation  into  atoms  are  given  "by  Herzberg(43).  Figures  for  both 
0°K  and  300°K  are  given  in  Table  2. 

(b)  Free  energies  and  equilibrium  constants.  The  free  energies  of  CO,  H*,  N»,  02 
and  graphite  have  been  interpolated  from  Wagman’s  data.  Those  of  0H(44),  110(45), 
H(46)f  n(47)  and  o(48)  have  been  corrected  to  new  values  of  the  fundamental 
constants(49)  using  similar  methods  to  those  used  by  Wagman.  For  C02  and  H20 
completely  new  calculations  have  been  made  by  the  present  author  using  the 
spectroscopic  data  given  by  Dennisan(50),  and  adding  the  correction  for  rotational 
stretching  term  given  by  Wilson(5l).  The  actual  expressions  used  for  C02  was 


*  =  16.0148  log  T  +  3.7640  +  +  .2385  T  .  10“5 

+  1.98718  { ln(  Qv),  +  ln(Qv)2  +  ln(Qv)s] 

+■  4.57565  log  (1  +  r) 

(8) 

where 

(qv),  =  2  exp  {-125j£ai  [n  -  .000223  (n2  +  n)]] 
n»o  a 

(9) 

(Q^)*  =  T  (n  +  1)  exp  {--^^  [n  -  .001945  (n2  +  n)]} 

(10) 

and 

(Qv)3  =  2  exp  {_2tl^ai  [n  -  .005264  (n2  +  n)]] 

nmO  A 

(11) 

the  upper  limits  for  n  being  chosen  such  that  the  exponentials  have  their  nrin-iimim 
value  there.  The  actual  values  are  of  no  importance  in  the  present  case  but  a 
finite  upper  limit  must  be  chosen  since  the  aims  to  infinity  diverge.  The  last 
term  1  +  r  is  composed  of  the  cross  terms  in  the  vibrational  levels,  terms  arising 
from  the  variation  of  moment  of  inertia  with  vibrational  quantum  numbers  and  also 
allows  for  the  existence  of  rotations  about  the  axis  of  symmetry  when  the  moment 
of  inertia  about  that  axis  is  made  non-zero  by  bending  vibrations.  The  complete 
expression  for  r  i3 

10*r  =  -  (1003  nj  +  2005  HJ) 

+  (1483  1ST  -  1150  H7  +  7850  5ST) 

♦  (2  n?  +  n|  +  62  n§) 


(12) 


(13) 

and  the  other  terms  can.  be  obtained  in  a  similar  manner.  To  calculate  the  sums 
of  these  series  the  (n2  +  n)  terms  in  the  exponentials  were  expanded  in  series  form, 


(3  nr  -  23  nr  ^  +  18  HS"  1ST) 

1 000 

+  (-  8199  Hi  na  +  31502  n,  n3  +  15823  xxg  n>) 

+  (34  n?  n |  +  496  n?  rtf  ♦  125  n|  n§) 

+  (258  n?  nj  H7  +  130  nT  n|  n7  -  498  nT  nT  n|) 

where 

n|  .  (q v)2  =  ^  n2(n  +  l)  exp  {-  ■9-6-1^26  [n  -  .001945  (n2  +  n)]] 
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so  giving  an  infinite  series  of  sums  of  the  farm 


N 

E  p(n)  sn 


n-o 


where  z  *  exp  {-  96l,26/Tj  eto.  The  error  in  taking  H  to  be  infinite 
negligible  and  so  each  sum  oould  be  written  down  as  a  polynomial  in  z  divided 
by  a  power  of  (l  -  z).  Ihe  work  was  made  much  easier  by  the  disoovery  of 
reourrenoe  formulae  far  successive  polynomials, 

Ihe  corresponding  expression  for  IfcO  is 

&  ■  -  8.1603  +  18.30268  log  i  ♦  *  2.33  T  .  1  <r* 


where 


and 


and 


♦  1.98718  [3n(Qr)|  ♦  3a(Qr)t  ♦  ln(Qv),l 
+  4.57565  log  (1  ♦  r) 

(Qv)t  *  E  exp  [-5252*^  [n  -  .01193  (n#  +  n)]] 
n«o  * 

(Qv)a  *  E*  exp  [-  [a  -  .01141  (nP  +  n)]] 

(Qr)»  *  #E^  exp  [-  [n  „  .01241  (n#  +  n)]] 


(14) 


(15) 

(16) 

(17) 


10*r  »  £1  (1000  +  72  np  -  126  5T  -  85  5T) 

-  (5019  nT  -  3101  5T  -  3345  ^) 

-  (133  hf  -  112  3*  -  104  3) 

-  (4 9  nT  aw  -  67  nT  -  34  n^  ^T) 


(18) 


These  expressions  are  somewhat  simpler  than  for  00t  but  the  shortness  of  Ihe  Q 
series  made  it  necessary  to  sum  (Qv),  term  by  term  for  T  >  4000  *K. 

Par  CH*  and  KBj»  new  calculations  were  made  using  the  moments  of  inertia  and 
vibrational  freqiencies  given  by  Dennison  in  a  rigid  rotator,  simple  harmonic 
oscillator  approximation  but  with  Wilson*  a  correction  terms  for  rotational 
stretching  added. 

Par  perfect  gases  the  formal  relation  for  thermodynamic  equilibrium  between 
one  group  of  gases  and  a  second  group  into  which  the  first  can  transform  by 
ohemioal  reaction  is  as  follows:-  Suppose  the  reaction  to  be 

lA  +  wB  +  nC  +  ..  jft  pX  +  qX  +  rZ  + 

where  l,a,n,  p,q,r  eto.  are  the  number  of  moles  of  eaoh  type  A,B,0,  X,X,Z  eto, 
being  lost  or  gained  in  the  reaction.  Let  the  partial  pressure  of  eaoh  type  in 
eqrilibrium  be  P(A)  eto, ,  then  the  ratio  of  the  products 


(p U)!1- |p(B)i“-  l?(o)!“* 

where  the  eqxilibcrium  constant  K  is  calculated  from  the  relation 

RT  In  X  *  {1  P*(A)  +  m  P°(B)  +  n  P*(C)  +  ...} 

-  [p  F*(X)  +  qP*(l)  +  r  P*(Z)  +  ...] 


(19) 


(20) 


where  the  free  energy  F°(a)  eto.  is  cal  ml  a  ted  for  one  mole  of  eaoh  substance  at 
temperature  T  but  at  unit  pressure.  Since  the  total  free  energy  of  one  mole  of 
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any  substance  contains  an  additive  term  of  R  times  the  natural  logarithm  of  its 
partial  pressure  we  see  that  our  equations  are  derived  directly  from  the  faot 
that  for  any  system  in  reversible  equilibrium  the  total  free  energy  is  at  a 
THirtnrtm,  in  the  tables  we  have  taioen  unit  pressure  to  be  one  atmosphere  and 
for  convenience  we  have  tabulated  (H?  -  F°  )/T,  where  (Ho  -  F°)  is  the  fall  in 
free  energy  of  one  mole  of  the  substance  when  its  temperature  is  raised  from  0°K 
to  T  at  constant  pressure  and  without  change  of  state.  If  we  assume  that  in  the 
reaoticn  we  cool  the  initial  substances  to  0°K,  carry  out  the  reaction  there  and 
then  raise  the  temperature  of  the  produo ts  back  to  T  then  we  can  use  the  faot 
that  at  0°K 

-  A  F?  =  -  A  H?  =  heat  of  reaction. 

We  can  therefore  calculate  the  free  energy  change  at  0°K  from  the  heats  of  formation 
at  that  temperature  given  in  table  2.  These  are  calculated  from  heats  of  formation 
at  room  temperatures  and  specifio  heats  over  the  range  down  to  0°K. 


As  an  example  oonsider 


£h2 


H 


for  which 


then 


-AH?  =  51620  cal. /mole.  Writing 

K,  =  p(H)/Vjp(Hk)] 


(o)  Heat  contents,  entropies  and  speoific  heats.  These  were  derived  from  the 
same  sources  as  were  the  free  energies,  the  data  for  OH  and  NO  being  correoted  in 
a  corresponding  maimer.  For  H,  N  and  0  it  is  sufficiently  accurate  to  take  the 
specifio  heat  to  be  that  of  an  ideal  monatomic  gas  namely  Cv  =  JB/2,  since  these 
gases  only  occur  in  very  small  amounts.  Data  for  CH*  and  NHs  were  calculated 
from  Dennison*  s  oons  tan  ts  as  before.  The  calculations  for  CO2  and  IfeO  were  based 
on  expressions  derived  from  those  for  the  free  energy  by  use  of  the  formulae 


H°  -  Hg 

T 


fa  -F» 


■] 


and 


C£  =  C  °  +  R  =  dH°/dT 


(21) 

(22) 


where  the  superscript  means  ideal  gas  conditions  i.e.  at  low  pressures.  The 
general  term  of  each  series  for  Qv  was  differentiated  and  in  differentiating  1  +  r 
use  was  made  of  the  fact  that  the  differential  coefficient  of  n*  can  be  expressed 
in  terms  of  n*-»i  etc.  The  resultant  expressions  are  too  lengthy  to  reproduce  here. 

The  entropy  at  one  atmosphere  pressure  was  calculated  in  each  case  from  the 
relation 


F  =  H  -  TS 


i.e.  3°  -  S^S.  .  (23) 

The  enthalpy  is  not  reproduced  but  I^T  can  be  obtained  if  required  by  subtracting 
entries  in  the  free  energy  table  from  corresponding  entries  in  the  entropy  table. 

It  was  considered  more  useful  to  tabulate  the  mean  specific  heat  at  constant  volume 


since  nearly  all  thermochemioal  calculations  oan  be  started  from  300°K. 

The  tables  have  been  entered  at  100°  intervals  so  that  linear  interpolation  is 
sufficient  in.  all  oases  except  perhaps  at  temperatures  below  2000°K.  The  logarithms 
of  all  the  equilibrium  constants  except  Ko  have  been  tabulated  since  these  logarithms 
have  almost  constant  first  differences.  Linear  interpolation  may  lead  to  an  error 
of  a  unit  or  two  in  the  third  decimal  place  -for  Ko  but  the  resultant  error  in  the 
calculated  temperature  will  amount  to  much  less  than  1°. 
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3.  Calculation  of  the  unoooled  temperature  of  explosion,  or  adiabatic  flame 

temperature. 

We  consider  all  solid  propellants,  including  colloidal  ones,  to  be  simple 
mixtures  i.e.  neglect  any  beat  of  mixing  there  may  be.  Liquid  propellants  are 
usually  not  mixed  before  firing  but  if  some  components  are  so  mired  e.g.  hydrogen 
peroxide  or  nitric  acid  diluted  with  water  the  heat  evolved  in  mixing  is  usually 
known  and  must  be  added  to  the  heat  of  formation.  We  oust  be  given,  or  we  must 
estimate,  the  fraction  ty  weight  of  each  ingredient  in  the  propellant,  multiplying 
each  fraction  by  the  beat  of  formation  of  unit  mass  of  that  ingredient  and  adding 
we  get  the  heat  of  formation  of  unit  mass  of  separated  ingredients,  to  which  must 
be  added  the  heat  of  mixing,  if  any,  to  obtain  the  heat  of  formation  of  unit  mass 
of  propellant.  $jr  similar  operations  on  the  other  columns  of  table  1  we  can 
obtain  the  number  of  gram  atoms  of  C,  H,  N  and  0  in  the  propellant  mixture. 

For  colloidal  propellants  and  monoprope llant s  composed  of  nitrated  organic 
liquids  suoh  as  nitrone  thane  or  ethyl  nitrate,  the  products  consist  almost  entirely 
of  the  five  gases  Ns,  COa ,  CO,  HtO  and  Ht,  the  proportion  of  inflammable  gases 
(CO  +  H*)  remaining  being  as  high  as  70  per  cent  in  some  oases.  With  bipropellant 
systems,  the  proportion  of  inflammable  gas  is  often  much  lower  and  when  hydrogen 
peroxide  is  used  as  a  malar  component  there  may  even  be  excess  oxygen.  Other 
bipropellant  systems  could  be  designed  to  produce  gases  containing  excess  oxygen, 
but  far  reasons  of  weight  these  are  never  used  in  practice  unless  the  oxygen  can 
be  picked  up  from  the  atmosphere  during  burning. 

Snail  amounts  of  inorganic  materials  are  often  included  in  colloidal 
propellants  and  large  amounts  in  plastic  propellants.  At  high  temperatures 
these  form  oxides  of  metals,  hydrides  of  halogens,  etc.  Sulphur  forms  H*S  plus 
a  small  proportion  of  SO. 

At  temperatures  below  1 ,000*0. ,  appreciable  amounts  of  methane  and  ammonia 
will  be  formed  in  cases  of  oxygen  deficit.  In  extreme  cases  solid  carbon  may 
separate  out.  At  temperatures  above  2,500°C.  appreciable  amounts  of  the 
dissociation  products  OH,  H  and  IK)  will  be  formed  and  at  still  higher  temperatures 
the  other  dissociation  products  0* ,  0,  N  and  C  will  begin  to  appear  in  that  order. 

m  the  very  short  time  interval  during  which  any  element  of  propellant  bums 
the  change  in  gas  pressure  will  be  very  small.  The  element  therefore  burns  under 
constant  pressure  and  makes  room  for  itself  by  doing  work  on  the  gases  previously 
formed.  For  explosion  under  constant  pressure  conditions  the  method  of  calculation 
of  the  temperature  of  suoh  an  element  is  as  follows.  We  make  a  rough  estimate  of 
what  that  temperature  will  be  and  then  calculate  the  composition  of  the  gases  in 
thermodynamic  equilibrium  at  the  estimated  temperature.  We  then  assume  that  the 
propellant  decomposes  completely  into  its  final  gaseous  products  at  constant 
pressure  and  constant  teisperature ,  the  hea-c  evolved  by  the  chemical  reactions 
being  stared  in  a  reservoir.  This  heat  is  then  taken  out  of  the  reservoir  and 
used  to  heat  the  gaseous  products  at  constant  pressure  up  to  the  final  temperature 
Symbolically:  - 

fT o 

-  A  %  (T  »  T*,  F  const.)  =  I  Cp  ( produo ts)  dT  (25) 

•'To 

where  Cp  is  the  specific  heat  of  the  gaseous  products  at  constant  pressure. 

Under  constant  volume  o auditions  each  element  as  it  is  produced  does  work 
not  on  the  external  world  but  on  the  other  elements  of  gas  produced  previously, 
which  will  be  ©Depressed  mare  or  less  adiabatioally.  The  first  element  to  be 

produced  will  suffer  qiite  a  high  subsequent  compression  and  so  may  be  expected 
to  reach  a  very  high  temperature  while  the  last  element  to  be  formed  will  remain 
at  temperature  Tg.  Convection  currents  or  turbulent  mixing  will  however  rapidly 
reduce  theee  temperature  differences  and  the  gas  will  take  up  a  temperature  T, 
given  by  the  fact  that  the  whole  of  the  heat  of  reaction,  at  constant  volume,  is 
transferred  to  internal  energy  of  the  gas  molecules,  i.e. 

fTl 

-AEr(T«T0,V  const.)  «  ,  Cy  (products)  dT 

% 
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Since  the  heat  of  reaction  at  oonstant  volume  usually  differs  try  only  a  few 
percent  at  roost  from  that  at  oonstant  pressure  we  have  roughly,  assuming  Tb/T© 
snwl  1 


Ti  ■  y  To 

where  y  is  a  mean  value  of  the  ratio  of  specific  heats  over  the  temperature 
range  up  to  T©  and  is  of  order  1,2  to  1.3. 

All  practical  methods  of  calculation  of  the  gas  ocooposition  proceed  by  a 
series  of  successive  approximations*  There  are  four  oases,  eaoh  requiring 
different  treatment.  The  first  three  are  for  calculations  at  temperatures 
below  2,700*0.  or  so  and  assume  that  the  amounts  of  dissociation  products  are 
relatively  small.  The  fourth  method  (d)  makes  no  suoh  assumption,  but  the 
total  pressure  cannot  be  predetermined.  The  first  three  oases  are  those  of 
(a)  oxygen  deficit,  (b)  coy  gen  excess,  and  (o)  oxygen  balance. 

Case  (a)  -  oxygen  deficit. 

We  start  with  1  gram  of  propellant  containing  (C),  (H),  (H)  and  (o)  gram 
atoms  of  eaoh  kind,  the  heat  of  formation  of  the  solid  propellant  being  h,  and 
we  wish  to  calculate  the  produots  formed  at  temperature  T0  under  a  total  pressure 
P0.  If  there  are  any  inorganio  elements  present  we  estimate  -the  produots  they 
form  first  and  subtract  the  necessary  amounts  from  our  gram  atoms  of  elements. 

We  then  assume  that  the  remainder  form  only  N*,  00,  CO*,  H*0  and  Expressing 

quantities  in  gram  moles  we  have 

(N*)  =  £(N) 

(CO)  =  (C)  -  (COa) 

(H*0)  *  (0)  -  (C)  -  (00a) 

(Ha)  =  £  (H)  -  (0)  +  (C)  ♦  (CO.) 

these  equations  expressing  the  equality  of  the  number  of  atoms  of  eaoh  t7P°  before 
and  after  the  reaction.  The  water  gas  equilibrium  gives  us: 

(CO)  x  (HkO)  =  Ko(00*)  x  (ife) 

this  equilibrium  being  independent  of  pressure.  Eliminating  (CO),  (H*0) 
by  means  of  the  previous  equations,  we  are  left  with  a  quadratic  equation 
(C02),  the  solution  of  which  is:- 

(CO.)  -  -  b  i  Vb*  ~  (0)  i(0)  -  (C))/|K,  -  li 

where  b  =  [(0)  +  Ko  (H)  -  (o)  +  (C)]]/2  (Ko  -  l) 

Only  one  of  the  two  solutions  far  (CO*  )  gives  positive  quantities  for  all 
produots.  We  can  now  write  down  the  number  of  gram  moles  of  each  of  the 
gases,  and,  by  proportion,  their  partial  pressures. 

The  next  stage  is  to  calculate  the  partial  pressures  of  the  dissociation 
products.  The  equilibrium  equations  should  preferably  express  each  dissociation 
produo t  in  terms  of  the  major  produots  only.  For  example,  we  should  use  the 
equation:  - 

n  P(0)  =  K,P  (Ife0)/p  (H,) 

v  " 

rather  than 

p(0)  =  (Ks/^Kl)  Mo*) 

since  p(0*)  will  not  be  known  very  precisely.  The  partial  pressures  of  the 
dissociation  products  oan  now  be  converted  into  gram  moles  using  tbs  same  faotor 
of  proportionality  as  for  the  major  gases.  These  dissociation  produots  are 
formed  by  reducing  the  amounts  of  the  major  gases  and  so  wa  must  subtract  the 
corresponding  amounts  from  the  (N) ,  (C),  (H)  and  (o)  available  for  the  major 
gases.  We  then  repeat  the  calculation  of  the  amounts  of  these  major  gases. 


and  (H*) 
for 

(27) 

(28) 

the 

major 
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Very  of  tan  ww  oan  guess  those  amounts  without  solving  •  the  quadratic  again;  a 
useful  rule  is  to  keep  the  amount  of  the  smallest  constituent,  usually  Ha,  the 
same  as  "before.  If  the  amounts  of  the  major  products  have  been  appreciably 
altered,  it  may  "be  neoessaxy  to  recalculate  the  dissociation  products  again. 

Having  calculated  the  gas  composition  at  the  given  total  pressure,  and 
estimated  temperature,  we  oan  calculate  the  heat  of  reaction  at  roam  temperature. 
This  is  simply  the  sum  of  the  heats  of  formation  of  the  gaseous  products  minus 
the  heat  of  formation  of  the  initial  propellant.  The  heat  required  to  raise  the 
gaseous  products  from  room  temperature  to  temperature  T0  oan  "be  calculated  using 
tables  of  mean  speoifio  heats  of  the  gaseous  products  over  this .  temperature  range. 
If  this  is  less  than  the  heat  available,  i.e.  the  heat  of  reaction,  the  true 
value  of  Iq  will  be  somewhat  higher,  and  vice  versa.  We  can  make  a  slight 
over-estimate  of  the  amount  we  have  to  raise  by  dividing  the  excess  heat 
available  by  the  speoifio  heat  at  temperatures  near  T0.  If  the  original  estimate 
of  T0  was  In  error  by  100*  or  more  it  may  be  neoessary  to  repeat  the  whole 
calculation  for  a  second  estimated  value  of  To  end  then  determine  the  true 
temperature  by  Interpolation. 


Case  (b)  -  oxygen  excess. 


This  oase  is  a  little  easier  than  (a)  because  we  have  no  quadratic  equation 
to  solve.  The  major  products  are:- 


and 


(N«)  =  ±(N) 

(CQ,)  =  (C) 

(Ife0)  =  i  (H) 

( O2 )  =  ±(0)  -(C)  -*(H) 


} 

] 


(29) 


The  minor  products  are  CO,  H2 ,  OH,  NO,  H,  0,  N  and  C,  usually  in  that  order 
of  magnitude.  The  equilibrium  equations  should  express  the  partial  pressures  of 
these  eight  dissociation  products  in  terms  of  the  four  major  gases.  Calculation 
then  proceeds  as  in  oase  (a). 

Case  (o)  -  oxygen  balance. 


In  this  case  we  have  only  three  major  products  and  the  calculation  of  the 
minor  products  is  more  difficult.  The  best  procedure  is  to  calculate  the  biggest 


minor  product  first,  namely  CO.  We  have 

(He)  ?  ^  {ior}  (co)  “  2^57  (co)  (30) 

an*  (Q«)  -  i f (CO)  ♦  (31) 

if  we  assume  that  CO  and  H2  are  the  only  minor  products.  We  also  have 

P(0s)  =  XSK«|(C02)/(C0)!2  (32) 

and  from  these  three  equations  vre  get,  since 

(CO.)  =  P„(0)/!(c)  +  i  (H)  +  £  (N)j ,  (33) 

(CO)*  =  W3£.(C)5!(C)  +  i  (H)  ♦  i  (N))/|2K.(C)  ♦  (H))  F0  (34) 


where  P0  is  the  total  pressure.  This  gives  us  a  first  estimate  of  (CO),  using 
which  we  oan  calculate  (H*)  from  the  exact  equation  and  then  (O2).  We  oan  then 
get  a  better  estimate  of  (CO)  by  inverting  equation  (32)  and  so  can  repeat  the 
cycle. 

Case  (d)  -  large  amounts  of  dissociation. 

In  such  cases  the  separation  into  major  and  minor  components  and  calculation 
by  successive  approximations  may  become  very  tedious  and  an  entirely  different 
method  of  approach  has  been  devised  by  Lutz! 52),  This  method  depends  on  the  faot 

that,  with  the  exception  of  NO,  the  partial  pressures  of  all  the  components  can  be 
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expressed  In  terms  of  two  parameters.  One  of  these  parameters,  namely  the  ratio 

*  *  (CO*)/[(CO)  +  (CO*)!.  (35) 

can  he  determined  from  the  initial  atomic  proportions  ty  use  of  a  universal  chart. 
The  other  parameter  is  the  partial  pressure  of  H<0.  A  value  met  be  given  to  this 


parameter  and  the  calculation  carried  right  through,  the  total  pressures  being 
obtained  as  the  sum  of  the  partial  pressures.  If  this  total  pressure  is  not  the 
value  required,  the  calculation  must  be  repeated  with  a  different  initial  value 
of  p(HzO).  The  method  is  therefore  one  of  trial  and  error.  Having  calculated 
the  gas  composition  for  the  required  temperature  and  pressure,  the  rest  of  the 
calculation  follows  normal  lines. 

An  outline  of  Iutz*s  method  is  as  follows: - 


We  write 


and 


a  =  x/(l  -  x)  =  p(COa)/p(CO) 
y  »  p(HgO) 

Pc  ■  P(CO*)  +  p(CO) 


PH  =  ?p(H.O)  +  2p(Ht )  +  P(OH)  +  p(H) 


?sr*^*K* 


Po  =  2p(C0a)  +  P(CO)  +  p(lfeO)  +  p(OH)  + 

+  Z&h)  +  P(X>)  +  P(l^) 

=  (1  +  x)  p0  +  4y-+  Zz  ^Koay  +  2  KlKga* 
+  ZsKoaj  +  p(NO^ 

=  (1  +  x)  Po  +  Pot  +  p(ND)  say 


(36) 

(37) 

(38) 


(39) 


(U>) 


If  we  divide  this  last  ecpaation  through  ty  Pg  we  havt 

E^=_Es  .  x  Si  „  £  +  igSl  (41) 

Since  p0,  pG,  eto. ,  are  proportional  to  (0),  (C),  etc.,  the  left-hand  side  is  a 
known  linear  function  of  x  for  any  given  composition.  If  we  omit  the  small 
correction  term  p(NO)  then  the  right-hand  side  is  a  function  of  x  (or  a)  only  for 
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a  given  value  of  y  for  a  given  temperature ,  since  the  Ks  are  functions  of  T  only. 

For  any  one  temperature,  therefore,  a  chart  can  he  prepared  on  'which  p</pg  is 
graphed  as  a  function  of  z  for  a  series  of  values  of  y.  The  curves  do  not 
differ  greatly  with  variation  in  y. 

For  any  given  mixture  the  straight  line  given  ty  the  left-hand  side  of 
equation  (41 )  oan  he  drawn  on  the  chart  and  where  this  intersects  the  p*/pjj  curve 
gives  the  correct  values  of  x  and  Pc/pg  for  a  given  y.  Given  the  values  of  x  and 
y,  we  oan  calculate  the  partial  pressures  of  all  the  hydrogen  compounds  and  so  pjj 
h y  emanation.  Ve  oan  then  calculate  pg  and  pjj  ty  proportion  and  so  the  partial 

pressures  of  the  remaining  components.  The  calculation  oan  he  shortened  slightly 
if  we  plot  a  set  of  curves  of  pg/y  on  the  same  chart.  These  again  are  functions 
of  x  and  y  only* 

We  oan  correct  approximately  for  2X)  formation  by  calculating  p(MO)/pg  and 
then  subtracting  this  from  the  left  hand  side  of  our  equation,  so  slightly  lowering 
the  line*  We  have 

{p(»D)A»KosJB  =  p(H«) 

*  **  "  *  ***>) 

giving  to  a  fair  approximation 

pCnoVph  -  AWWtoJ 

A  second  approximation  is  not  usually  necessary. 

The  method  requires  the  preparation  of  a  chart  for  each  value  of  T,  say  every 
100°,  over  the  range  of  interest  hut  these  charts  oan  he  prepared  once  and  for  all* 

% 

Fig*  1  shows  a  simplified  diagram  for  T  =  3*000*K.  It  oan  he  seen  that  the 
curves  do  not  differ  greatly  for  quite  large  changes  in  y,  so  that  we  do  not  have 
to  (press  the  initial  value  of  the  partial  pressure  of  H*0  very  exactly.  Of  oouras, 
in  practice  curves  would  he  drawn  an  the  diagram  for  a  considerable  number  of  values 
of  y.  Far  temperatures  below  2,500°K  the  curves  coincide  for  all  practical  values 
of  y  and  x  and  the  method  then  degenerates  into  a  graphical  method  of  solving  the  . 
quadratic  equation  (27)  for  (CO*). 

Having  calculated  the  gas  composition  at  the  given  total  pressure  or  volume 
and  the  assumed  temperature  we  use  the  table  of  mean  specific  heats  to  calculate 
the  heat  required  to  raise  that  gas  mixture  to  that  temperature  from  300*K.  This 
should  equal  the  heat  of  reaction  available  for  the  purpose,  namely  the  heat  o if 
formation  of  the  gaseous  products  which  oan  he  calculated  from  table  2,  less  the 
heat  of  formation  of  the  initial  propellant,  the  whole  calculation  being  carried 
out  for  constant  pressure  or  constant  volume  as  the  ease  may  be.  If  equality  is 
not  attained  the  whole  calculation  must  be  repeated  for  another  assumed  temperature 
T0  or  Ti .  For  an. experienced  operator  two  such  calculations  are  visually  sufficient. 

Finally  mention  should  be  made  of  a  very  simple  approximate  method  of 
calculation  developed  by  Hirsohfelder  and  Shermam  53 ) .  This  oan  only  be  used  as  it 

stands  for  mixtures  deficient  in  OJQrgen  and  for  vary  small  amounts  of  dissociation, 
so  that  it  is  mainly  of  value  for  gun  propellants.  For  values  of  Ti  between  2000 
and  3000°K  the  three  equations  to  be  used  are 

n  »  8  3^4  (42) 

Oy.  ■  2  *i  Cvi  (*3) 

and  T,  =  2500  +  (L  x±  Ii)/CT  °K  (44) 

where  x^  is  the  fraction  by  weight  of  the  i  th  constituent  of  the  propellant  and 
n^,  Cyi  and  are  tabulated  for  all  the  usual  propellant  ingredients.  For 
temperatures  over  3,000°K  dissociation  reduces  Ti  appreciably  and  a  more  elaborate 
formula  is  given  for  Ti  in  such  oases.  Strictly  speaking  their  table  should  be 
modified  slightly  to  bring  it  into  line,  with  the  more  accurate  thermoohemiqal  data 
now  available  but  the  approximate  nature  of  their  method  renders  this  unnecessary. 
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4.  Corrections  for  gas  imperfections. 

Under  rocket  conditions  i.e.  temperatures  greater  than  1 ,000°C  and  pressures 
less  than  100  a  too  spheres ,  wb  may  treat  the  gaseous  products  as  a  mixture  of 
perfect  gases.  For  gun  conditions  however  the  pressures  are  of  such  a  magnitude 
that  fairly  simple  corrections  are  both  necessary  and  sufficient.  Gun  calculations 
invariably  come  under  case  (a). 

•  These  pressure  corrections  are  treated  fully  by  Comer (3)  tut  for  convenience 
a  short  abstract  of  his  paper  is  given  here.  in  the  first  place  the  equation  of 
state  of  such  a  gas  mixture  is  represented  with  sufficient  aocuracy  by 

F7  =  nRT  (l  +  BA  +  nC/V2)  (45) 

where  B  and  C  are  additive  functions  of  the  constituents  of  the  gas  mixture,  e.g. 

B  =  (COa)  •  B(C02)  +  (CO)  .  B(C0)  +  ...  (46) 

and  similarly  for  C.  The  water  gas  equilibrium  constant  is  multiplied  by  the 
factor 

exp  (-  n  ABA  "  n2  AC/2V* ) 

where  AB  =  B(C0)  +  B(H20)  -  B(C00  -  B(Ha)  (47) 

with  a  similar  expression  for  AC.  The  only  other  equilibrium  constants  that  need 
correcting  are  Ke  and  K3  and  it  is  sufficient  to  multiply  each  by  exp  (-20  nA)» 

The  internal  energy  at  some  temperature  T  ua n  be  written 

E  -  E300  =  S(C02)  Eo(C02)  +  y  E(C02)  •  B,(C0a)  +  |^y  2(C02)  E2(C02)  (48) 

where  -the  summation  signs  indicate  sums  over  all  the  constituents  of  the  gaseous 
mixture.  The  first  term  is  the  energy  of  the  gases  assumed  perfect  relative  to 
the  same  perfect  gases  at  300°K,  and  the  correction  terms  Ei  and  E2  for  the  various 
gases  are  given  by  Corner.  The  nett  correction  to  E  is  always  negative.  It 
should  be  noted  that  no  correction  is  to  be  applied  to  the  heat  of  reaction,  which 
has  been  corrected  to  perfect  gas  conditions  in  preparing  our  tables  of  heats  of 
formation. 

For  imperfect  gases  it  is  best  to  work  in  terms  of  partial  densities,  not 
partial  pressures,  the  equilibrium  equations  are  then  used  in  the  form 


(H)  =  K,  ^i(H2)  V/82.06  T]  (50) 

(OH)  =  Kb  (HaO)  ^[v/82.06  T^)]  exp  (-20  i\A}  (51 ) 

<52> 

(N)  =  K4  V{(lfe)  V/82.06  T}  (53) 

(°)  =  £5  82.66  T 

(O’)  (55) 


where  the  bracketed  quantities  are  moles  per  gram  of  mixture,  V  is  the  given  volume 
of  one  gram  of  mixture  and  n  the  total  number  of  moles  of  gas  in  one  gram  of  mixture 
For  use  in  these  equations  it  is  usually  sufficiently  accurate  to  calculate  n 
assuming  no  dissociation,  i.e.  to  put 

n  =  (C)  +  -5-  (H)  +  y  (N)  (56) 
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where  in  this  equation  the  quantities  on  the  right  hand  side  are  the  total  number 
of  gram  atoms  of  C,  H  and  N  per  gram  of.  propellant. 

Having  completed  the  calculation  of  T<  we  can  compute  B  and  C  in  the  equation 
of  state  and  so  the  pressure  P.  But  in  internal  "ballistic  calculations  for  gun« 
it  is  far  more  convenient  to  use  an  eqiation  of  state  of  the  farm 

P(y  -  Tj)  =  hRT  (57)  * 

where  the  co-volume  tj  is  taken  to  be  constant.  If  we  choose  a  suitable  value  for 
V  near  the  minimum  value  that  occurs  in  the  gun  then  we  calculate  a  value  of  77 
that  gives  the  correct  pressure  at  that  gas  density  and  very  nearly  the  correct 
value  at  all  densities  since  at  low  pressures  the  effect  of  corrections  for  gas 
imperfection  diminish.  It  is  customary  to  take  V  «  5  0.0./ gjn.  and  so  we  have 
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25B  +  5nC 


(58) 


Since  the  C  term  is  relatively  small  and  B  -  1 ,  we  find  n  to  be  of  the  order  of 
80#  of  B. 

Corner  also  gives  a  rough  formula  for  estimating  77,  namely 

77  =  1.18  +  6.9  (C)  -  11.5  (0)  (59) 

where  the  quantities  on  the  right  hand  side  are  gram  atoms  per  gram  of  propellant. 

For  calculating  the  effective  value  y  of  the  ratio  of  specific  heats  for  use 
in  gun  calculations  we  must  allow  for  heat  liberated  by  chemical  reactions  as  the 
gases  cool,  assuming  thermodynamic  equilibrium  to  be  maintained  as  is  almost 
certainly  the  case  at  the  high  temperatures  and  pressures  obtaining  in  the  gun. 

For  this  purpose  we  assume  that  when  the  base  of  the  projectile  reaches  the  nuzzle 
T  has  fallen  to  0.6  times  its  uncooled  value,  and  V  haw  risen  to  10  c.o./ga.  We 
then  calculate  the  gas  composition  in  equilibrium  under  those  conditions,  an  eaqjr 
calculation  since  there  will  be  no  dissociation.  The  difference  between  the  heat 
of  reaction  available  and  the  heat  content  of  the  gases,  relative  to  ideal  conditions 
at  300°K.  gives  the  energy  lost  to  external  work.  If  m  call  this  energy  loss  AE 
then  we  have 


y  =  1t0^I  (60) 

where  y  is  the  effective  mean  value  that  we  require. 

5.  The  use  of  entropy  tables  in  calculating  rocket  thrusts. 

m  a  gun  we  only  get  adiabatic  conditions  after  the  charge  is  all  burnt  and 
it  is  sufficient  to  use  the  mean  value  of  y  as  defined  above,  but  as  a  general  rule 
the  use  of  a  constant  y  is  not  suitable  for  calculating  the  thrust  obtainable  from 
a  rooket. 

The  total  thrust  given  by  a  rocket  motor  is  the  resultant  farce  of  the  gas 
pres  sura  acting  over  all  surfaces  of  the  motor.  If  Pe ,  the  pressure  at  the  exit 
of  the  nozzle,  equals  Pa.  the  atmospheric  pressure,  then  the  atmosphere  produces 
no  resultant  thrust  on  the  motor.  Otherwise  there  is  an  additional  thrust  equal 
to  (Pe  -  Pa)  multiplied  by  the  nozzle  exit  area.  The  main  component  of  thrust  is 
equal  to  the  rate  of  generation  of  momentum  i.e.  the  product  of  the  gas  exit 
velooity  multiplied  by  the  mass  discharged  per  second. 

The  calculation  of  the  pressure-temperature  relation  in  the  adiabatic  expansion 
through  the  nozzle  is  more  difficult  if  there  are  appreciable  quantities  of 
dissociated  gases  in  the  chamber  since  these  may  recombine  as  their  temperature 
falls  while  they  go  through  the  nozzle.  Thera  is  some  doubt  about  the  speed  of 
this  recombination,  but  at  the  very  high  temperatures  at  which  an  appreciable  amount 
of  dissociation  occurs  chemical  reactions  are  extremely  fast  and  it  seems  reasonable 
to  suppose  that  chemical  equilibrium  is  maintained  at  least  as  far  as  the  nozzle. 
This  means  that  most  of  the  recombination  will  take  place  and  so  we  shall  not  be 
far  wrong  if  we  assume  it  to  be  complete. 
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Ill  calculating  the  pressure -temperature  relation  there  is  some  advantage  to 
"be  gained  if  we  make  use  of  the  condition  of  constant  entropy.  We  must  first 
calculate  the  gas  composition  at  the  temperature  Tc  in  the  chamber  and  at  a 
total  pressure  P0.  If  S°(C(V)  is  the  entropy  of  one  mole  of  CQ*  at  the  given 
temperature  and  a  partial  pressure  of  one  atmosphere  then  the  entropy  per  gram 
of  our  gas  mixture  at  temperature  T  and  pressure  P  is 

S  =  2  (CO*)  S°(CO*)  -RE  (00*)  P(CO,)  (6l) 

where  P(CO*)  is  the  partial  pressure  of  CO*.  Taking  the  gases  to  he  perfect  we 
have 


P(COa)  =  P  .  (CO*)/n 

and  so 

S=S°  -  Rn  In  P  -  Rn  £  -^2*1  in  lS2tl 

n  n 


(62) 


where  we  have  written  S°  for  the  first  sum. 

Usually  the  final  pressure  Pe  is  stipulated  and  so  we  must  calculate  the  gas 
composition  and  entropy  for  at  least  two  temperatures  at  that  pressure.  The 
correct  temperature  can  then  he  interpolated. 

If  no  changes  in  chemical  composition  occur  then  the  last  term  in  equation 
(62)  remains  constant  and,  since  we  are  only  interested  in  differences  of  entropy, 
it  need  not  he  calculated.  For  an  adiabatic  expansion  we  then  have  6S  =  0 


i.e.  Rn  6  (log  P)  =  6  S°  =  6  T  (63) 

where  n  Cp  =  2  (CO*)  Cp(CO*)  (&».) 

We  can  calculate  the  increase  in  kinetic  energy  of  the  gases  from  the  relation:  - 

ib  (v2)  =  -nCp  6T  (65) 

or  i  6  (v2)  =  RTn  6  (In  P)  (66) 


If  chemical  changes  occur  during  the  expansion  then  for  no  change  in  total  entropy 
we  have 


i  b  (v2)  =  -  n  6  H  (67) 

where  we  have  generalised  -  6  H  to  include  not  only  -  Cp  6T,  the  purely  thermal 
energy  liberated,  hut  also  the  heat  liberated  by  the  chemical  reactions  occurring 
during  the  small  temperature  drop.  Some  writers  use  a  generalised  H  which 
satisfies  equation  (67)  and  which  is  defined  by  the  relation;  - 

nH  =  2  (CO*)  H'(CO*)  (68) 

r T 

where  H'(CO*)  =  j  Cp(CO*)dT  -  (H.F.  of  CO*  at  0°K)  (69) 

•'0 

H*  will  then  have  normal  values  for  elementary  gases  such  as  H*  and  N*  but  largo 
negative  values  for  such  gases  as  CO* ,  CO  and  H20.  With  these  generalised 
enthalpies  equation  (67)  will  be  true  whether  chemical  changes  take  place  or  not. 

For  any  given  gas  mixture,  diagrams  connecting  total  entropy,  generalised 
enthalpy,  pressure  and  temperature  can  be  constructed  in  precisely  the  same  way 
as  for  steam.  Such  diagrams  are,  however,  only  accurate  for  tliat  particular 
mixture  and  so  are  rarely  worth  constructing.  For  a  series  of  mixtures  a  series 

of  diagrams  would  have  to  be  constructed  but  if  the  whole  series  had  eqial  values 
of  some  parameter,  such  as  the  molar  fraction  of  tria tonic  gases  in  the  mixture, 
then  the  diagrams  would  be  nearly  the  same  and  for  rough  calculations  a  general 
diagram  could  be  used.  A  good  deal,  of  calculation  on  these  lines  has  been 
carried  out  by  German  scientists. 
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Equation  (66)  is  always  true  and  this  equation  does  not  contain  Cp.  Hence 
for  a  small  pressure  drop  the  gain  in  kinetic  energy  of  the  gases  is  independent 
of  whether  or  no  the  heat  liberated  is  inoreased  try  chemical  reactions ,  provided 
we  start  from  the  same  initial  conditions.  Since  most  of  the  chemical  reaction, 
if  it  does  occur,  will  occur  during  the  first  200°  fall  in  tenperature  its  effeot 
on  the  thrust  will  be  quite  small.  Unless,  therefore,  the  dissooiation  is 
Riffioient  to  lower  the  chamber  temperature  try  more  than  100°  or  so  the  thrust 
can  be  calculated  with  sufficient  a o curacy  for  most  purposes  by  a  simplified 
calculation  assuming  that  no  chemical  reaotions  take  plaoe  in  the  nozzle. 

Once  we  have  calculated  the  temperature  and  composition  of  the  gases  under 
a  given  pressure  PQ  in  -the  combustion  chamber,  by  one  of  the  methods  outlined  in 
seotion  3,  the  simplified  calculation  for  the  adlabatio  expansion  is  very  eaey. 

In  an  adiabatic  expansion  the  total  entropy  S  remains  oonstant  but,  since  we  are 
assuming  the  gas  ccoposition  to  remain  oonstant,  the  third  term  on  the  right-hand 
side  of  equation  (62},  whioh  is  in  faot  the  entropy  of  mixing  of  the  gases,  remains 
oonstant  and  so  need  not  be  calculated.  Under  these  conditions  in  an  adlabatio 
expansion  3°  -  Rn  log  P  remains  oonstant.  S°  is  a  function  of  temperature  only 
and  can  be  readily  calculated  from  tables  for  a  given  gas  mixture.  The  procedure 
is  to  calculate  S6  at  the  combustion  ohamber  temperature  T0  and  than  at  soma  lower 
temperature  Tf •  Since  S  remains  oonstant  we  have:- 

S°(T0)  -  S°(Tk)  =  Rn  (In  Pc  -  In  P,)  (70) 

from  which  we  can  calculate  the  pressure  P*  of  the  gas  corresponding  to  the 
temperature  T2.  By  choosing  two  or  more  values  of  T2  and  then  interpolating 
to  can  calculate  the  value  of  Te  corresponding  to  the  nozzle  exit  pressure  PQ. 

The  kinetic  energy  of  the  gases  issuing  frcm  the  nozzle  will  equal  the 
deorease  in  their  enthalpy  H,  which  is  a  function  of  temperature  only,  for  a 
given  gas  mixture.  H(T0)  has  been  calculated  already  in  evaluating  T0.  We 
then  calculate  H(Te)  and  to  have: - 

*vf  =  H(To)  -H(Xe)  (71) 

The  speoifio  impulse  is  the  gas  momentum  produced  per  unit  mass  of  propellant 
consumed  and  if  the  products  are  entirely  gaseous  this  is  simply  ve  or  Vg/g  in 
practical  units.  The  thrust  on  the  motor  is  this  quantity  multiplied  by  the  mass 
rate  of  consumption  of  propellant.  We  are,  of  course  assuming  expansion  to 
atmospheric  pressure  in  the  nozzle.  The  specif io  impulse  increases  slowly  with 
chamber  pressure  and  striotly  speaking  should  be  corrected  to  some  standard 
chamber  pressure,  say  20  atmospheres.  Since  in  practice  the  range  of  ohamber 
pressures  used  is  not  very  wide  this  correction  is  usually  emitted. 

6.  Numerical  examples. 

To  illustrate  the  various  methods  of  calculation  to  give  two  examples,  one 
for  a  gun  propellant  with  very  little  dissociation  of  the  gaseous  products  and 
one  for  a  highly  dissociated  rocket  propellant. 

(a)  Calculation  for  S.C.  in  a  oonstant  volume  of  5  c.c./gn« 

S.C.  may  contain  seme  .2$  by  weight  of  water  and  say  ,5%  of  inert  material, 
such  as  chalk,  iron  oxide,  silica  etc.  which  occur  as  impurities  in  the  nitro¬ 
cellulose.  Taking  nominal  proportions  of  the  dry  ingredients  we  get  the 
following  composition  per  gm. 

NO  .48902  gms. 

NO  .41417 
Carb.  ,08982 
IfeO  .00200 
Inert  .00499 

First  of  using  Hirsohf elder ' s  simple  method  to  have,  assuming  a  speoifio 
heat  of  .3  oal./ga.  °K  for  the  inert  material  which  absorbs  heat  but  contributes 
nothing  to  the  pressure, 


of  12.3JS  N  (actual) 

Gal.val.  (calculated)  934  oal./gn. 
H.F.  477  cal./ga. 


1 
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n  =  .04233  gm.  molee/gm. 

Cy  =  .35013,  y  =  1  +  1.987  ry^  =  1.240 
Ti  *  2500  +  21 7.55/.35013  =  3121  °K. 

The  temperature  will  be  a  little  lower  on  account  of  dissociation,  say  3100°K, 
which  we  shall  use  as  our  first  estimate  far  the  method  of  case  (a).  The  atomic 
composition  is 

(C)  =  .022107  gm.  atonm/gm. 

(H)  =  .02998^ 

(N)  =  .010435  n  ~  (C)  +  £  (H)  +  £  (N)  =  .042317 

(O)  =  .034589 

[Check  .00499  +  12.010  x  .022107  +  1.008  x  .029984  +  eto.  «  .999996  @ns.] 

(CO)  =  .022107  -  (C02) 

(H*0)  =  .012462 -  (CO*) 

(H*)  •=  .002530  +  (CO*) 

and  K  =  7.923  exp  [.2674  +  .OI83]  s  10.543 

Hence  (CO,)  =  -  .003209  +  /. 0000391 655  =  .003049 


Our  calculated  products  are 

*  now 

(CO*)  =  .003049  gm. 

moles 

(H2O)  =  .009413  - 

•• 

I” Check  (ifeO)  •  (CO)  =  .00017939 

(co)  =  .019058  .. 

M 

L  K  (Ha)  .  (CO, )  =  .00017934 

(Ik)  =  .005579  - 

•• 

(n2)  =  .005218  .. 

— 

Inert  =  .00499  ©ns. 

and  from  these  we  calculate  the  following  byproducts,  using  equations  (50)  to  (52) 

(H)  =  .2121  V[ .005579/(1 6.41 2  x  3100)]  =  .000070 

(OH)  =  .07695  x  .009413/1.1845  ^16.412  x  3100  x  .005579] 

=  .OOOO36 

and  (HO)  =  .008519x  .009413  V].0052l8]/l.l845  x  .005579  -/[16.412 x3100] 

=  .0000C4 

the  other  byproducts  being  negligible.  Subtracting  these  from  the  gm.  atoms 
available ' for  the  main  products  we  have 

•  (C)  =  .022107  gm.  atoms/ gm. 

(H)  =  .029878  ..  ..  ..  n  =  (C)  +  i  (H)  +  £  (N)  +  .000110 

(N)  ‘  =  -.010431  -  .042372 

(O)  =  .034529  ~ 

(CO)  =  .022107  -  (CO2) 

(HaO)  =  .012422  -  (CO2) 

(H*)  =  .002517  +  (CO2) 

and  K  =  7.923  exp  [.2678  +  .0183]  =  10.547 

We  can  now  find  the  major"  products  by  solving  the  quadratic  or  by  trial. 

.The  changes  are  so  small  .as  to  render  recalculation  of  the  major  products 
unnecessary  and  we  have  for  the  complete  composition  of  1  gm. 
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Subst. 

Amount 

»  m  ■  <TOT-i 

H.P. 

4 

(CO*) 

•003046  gm.  moles. 

94052  cal  ./mole. 

11.562  cal  ./mole.  °K 

(h*o) 

.009376  .. 

57504  - 

— 

9.288  .. 

•• 

(CO) 

.019061  ..  .  .. 

26711  - 

— 

6.312  .. 

•• 

(Ha) 

.OO5563  - 

- 

5.881  .. 

— 

(N.) 

.005216  « 

- 

6.24-5  •• 

•• 

(H) 

,000070  ». 

-51794  - 

— 

2.981  .. 

M 

(OH) 

.OOOO36  ». 

-  9903  - 

— 

5.988  .. 

— 

(HO) 

.000004  *.  «• 

-21600  .. 

— 

6.445  - 

— 

Eiert 

.00499  gm. 

«■ 

0.3  oal./gn.  °K 

n  s 

.042372  ,  H.P. 

-  1330.7  oal 

•/gpu 

,  C|  =  .30979  oal./gn.  °K 

The  nett  heat  available  from  the  reaoticn  is  1331  -  477  =  854  cal./gn. ,  the 
reduction  from  the  cal.val.  being  due  to  the  'water  now  being  gaseous.  The  heat 
required  to  heat  the  produots  at  constant  volume  but  at  very  low  pressure  from 
300°K  to  31 00°K  is  Eo  =  .30979  x  2800  =  867  cal./ gm.  but  to  this  must  be  added 
the  negative  terms  due  to  gas  imperfections  as  explained  in  section  4.  Using 
Corner* s  tables  vre  have 

Sjioo  “  Esoo  =  867  -  .0084m  x  428  +  (.0084744)2  x  18404- 
=  867  -  3.6  +  1 .3  =  865  cal./gn. 

This  exceeds  the  heat  available  by  11  oal./gn.  Using  table  7  we  find  that  the 
specific  heat  of  the  products  at  3l00°K  is  .356.  The  true  value  of  Ti  will 
therefore  be  a  little  higher  than  3100  -  11/.356  =  3069°K  because  the  heat  of 
reaction  will  be  a  little  greater  at  that  temperature.  The  main  difference  is 
that  the  byproducts  are  reduced  by  a  factor  of  about  0.7  for  a  drop  in  temperature 
of  1 00°K  and  the  heat  absorbed  in  dissociation  at  3100°K  is  6.4  cal  ./go.  in  the 
present  case.  A  better  estimate  of  Tt  is  therefore  3100  -  10/. 356  =  3072°K. 

At  that  temperature  we  have  by  interpolation  of  table  4,  Ki  =  .1952, 

Kg  =  .06931 ,  Kj  =  .007524  and  since  the  changes  in  the  major  products  will  be 
small,  we  have  for  the  new  values  of  the  byproducts 

(H)  =  .000070  x  .195^.2121  =  .000065  gm.  moles, 
similarly  (OH)  =  .000033,  (N)  =  .000004-  gm.  moles. 

Allowing  for  the  increase  in  the  major  products  of  -jr  [.000005  +  . 000003 ]  and  a 
decrease  in  the  minor  products  of  twice  that  amount  we  have  n  =  .04^368  gm./ moles/ 
gn.  Interpolating  K<>  we  have 

K  =  7.868  exp  [.2678  +  .0184]  =  10.4-75 
and  (H*0)  =  .012425  -  (CO*) 

(CO)  =  .022107  -  (C02) 

(Hg)  =  .002518  ♦  (C02) 


We  can  find  the  major  products 
is  as  follows. 

1  by  trial  and  the 

final  composition  at  3072°K 

Subst. 

Amount 

H.P. 

(CO*) 

.OO3056 

gn. 

moles. 

CM 

oal./mole. 

11.549 

cal  ./mole  °K 

(H*0) 

.009369 

•• 

•• 

57504 

—  — 

9.266 

•• 

— 

(CO) 

.019051 

•• 

•• 

26711 

—  — 

6.306 

— 

M 

(H*) 

.005574 

— 

— 

- 

x 

5.871 

— 

— 

(N*) 

.005216 

— 

- 

6.238 

— 

— 

(H) 

.OOOO65 

— 

— 

-51794 

—  — 

2.981 

— 

— 

(OH) 

.OOOO33 

— 

•• 

-  9903 

—  w 

5.979 

— 

— 

(NO) 

.000004. 

— 

•• 

—21 600 

—  — 

6.439 

— 

— 

Inert 

.00499  gm. 

- 

.3 

cal./gn.  °K 

n  =  .042368 

$ 

H.P. 

*  1331. 

,3  cal./gn. 

,  C®  =  .30942  cal./gn.  °K 
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Henoe  Ejot*  “  Eloo  =  855  oal./gm.  since  the  pressure  correction  will  be  -2  oal./gm. 
as  before*  The  heat  available  is  854  cal./gm.  so  the  true  temperature  is  3070®K, 
the  composition  being  unchanged  from  the  above.  The  force  constant  is 

X  =  *5384  n  T,  =  70.05  tons/sq.in./gm./o.o. 
or  n  1938  inch  tone/lb. 

The  corresponding  figure  calculated  by  the  simple  method  is  higher. 

To  calculate  y  we  must  calculate  the  gas  composition  at  0.6  x  3070  =  1842  #K 
and  at  a  volume  of  10  o.o./gm.  There  will  be  no  byproducts  and  the  water  gas 
equilibrium  constant  is 


K  =  4.177  exp  {.1413  +  .0076]  «  4.848 

The  products  are  now 

Subst. 

Amount 

H.F.  c| 

(CO*) 

.OO4337  gm.  moles. 

94052  cal  ./mole.  10.688 

cal./ mole.  #K 

(H*0) 

.008125  - 

57504  -  -  8.042 

••  — 

(CO) 

.01 7770  —  w 

26711  -  «  5.916 

H  M 

(Ha) 

.006867  —  — 

5.375 

—  M 

(Na) 

.005218  « 

5.833 

M  M 

Inert 

.00499  fiP*. 

.25 

oal./gm.  °K 

n  =  .042317  ,  H.F.  a  1350  cal./gm.  ,  ■  .28542  cal./gm.  °K 

Et  *4*  -  IS 00 

n  440.1  -  4.9  +  .4  =  436  oal./gm. 

Heat  available 

=  1350  -  477  =  873  cal./gm. 

Hence  the  nett  loss  of  internal  energy  when  the  temperature  falls  from  3070  to 
1842°K  is  437  cal./gm.,  so  over  that  range  Cy  =  .3559,  n  =  .042343  and  so 
y  =  1 .2364.  This  compares  with  1 .240  as  calculated  by  the  simple  method. 

The  second  and  third  virial  coefficients  in  the  equation  of  state  calculated 
at  30700 K  are 

B  =  1.1170  c.o./gn. 
n  C  =  .2570  (c.o./gm.)a 

Hence  from  equation  (58)  i)  =  .947  o.o./gn.  Allowing  for  the  volume  of  the  inert 
solids  would  inorease  this  by  about  two  units  in  the  third  decimal  place,  we  may 
therefore  round  off  this  figure  to  .950  c.c./@n. 

The  value  of  tj  calculated  for  V  =  10  c.c./gm. ,  T  =  1842°K  is  1.002  o.c./gm. 
but  since  the  co -volume  is  much  less  important  here  it  is  sufficient  to  take  the 
constant  value  of  .950  c.c./gm.  The  simple  formula  equation  (59)  suggested  by 
Comer  gives  rj  =  .935  c.c./gm.  which  is  sufficiently  accurate  for  most  purposes. 

(b)  Calculation  for  a  1:3  mixture  by  weight  of  decane  and  liquid  oxygen  under 

a  constant  pressure  of  50  atmospheres. 

Decane  has  been  ohosen  as  a  representative  paraffin.  The  mixture  gives 
values  of  Tq  of  the  order  of  3600°K  and  so  large  amounts  of  dissociation.  We 
shall  therefore  start  by  vising  Lutz's  graphical  method.  The  atomic  composition 
of  this  mixture  is 

(C)  =  .017572  ©a.  atomVgm. 

(H)  =  .038657 
(0)  =  .046875  •« 

H.F.  =  199  cal./gm.  at  C.P. 
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We  shall  use  Lutz'B  method(*“7  to  find  the  gas  composition  at  3600°K.  Sinoe 
p(HzO)  must  "be  of  order  20  atmospheres  we  see  from  his  graph  that  x  nust  lie  between 
.4  and  .5,  but  since  his  equilibrium  constants  differ  from  ours,  particularly  far 
OH,  we  shall  draw  our  own  curves  of  Po/pg  for  y  a  20  and  y  =  l6  atmos.  and  for 
values  of  x  extending  only  from  .3^  to  .50*  Putting  in  our  own  equilibrium 
constants  we  have 

Pet  (y  **  20)  =  20  +  4*808  /a  +  1,410  a  +  9.102  a* 

Pa  (y  =  16)  =  16  +  4*300  Va  1*410  a  +  9*102  a* 

Pjj  (y  ®  20)  =  40  +  4*808  Va  +  l.llJ/Va  +  4.575/a 

Ph  (y  =  16)  =  32  +  4.300  Va  +  *9W/Va  ♦  4.300/a 

where  a  =  x/0  “  x).  The  intersecting  line  is 


.75803 

-  .45456  x 

and  the  points  of  intersection  give 

✓ 

y  a  20,  x  =  .4135, 

PH  =  51.85, 

a  a  ,7050 

y  a  16,  X  =  .3977, 

PH  =  42.26, 

a  a  .6603 

Prom  cur  equilibrium  constants  we  can  now  calculate  two  trial  solutions,  the 
quantities  in  ga.  moles  being  calculated  from  the  partial  pressures 


Gas 

p  atmos. 

go.  moles/gm. 

p  atmos. 

gn.  molee/gn. 

H^0 

20 

.014912 

16 

.014636 

H. 

3.244 

.002419 

2.771 

.002535 

OH 

4.037 

.003010 

3.495 

.003197 

H 

1.324 

.000987 

1.224 

.001120 

02 

2.262 

.001686 

1.984 

.001815 

0 

•99k 

.000741 

.931 

.000852 

CO* 

.007267 

7.640 

.006989 

CO 

C.823 

.010305 

■lliflg 

.010583 

Total 

55.430 

.041327 

45.615 

.041726 

C 

23.569 

.017572 

19.210 

.017572 

H 

51.849 

.038657 

42.261 

.038658  ‘ 

0 

.046874 

.046876 

j  Lit  *3 

By  repetition  we  can  find  the  correct  value  of  y  which  turns  out  to  be 
17*78  atmospheres.  By  simple  linear  interpolation  of  the  two  trial  solutions 
we  find  y  =  17*79  atmospheres  and  linear  interpolation  of  the  two  gas  compositions 
as  given  in  gm.  moles/ gm*  leads  to  quantities  in  error  by  about  3  in  the  last 
decimal  place.  These  errors  combine  to  give  an  error  of  1  cal./ga.  in  the  heat 
of’ reaction,  a  quite  negligible  amount,  so  that  there  is  in  practice  no  need  to  • 
calculate  another  solution  provided  one  value  we  have  chosen  for  y  is  not  mare 
than  10$  in  error;  This  limit  applies  when  Tq  =  3600°K  and  considerably  wider 
limits  can  be  allowed  at  lower  temperatures. 

The  exact  solution  is  as  follows 


/ 
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Gas  gm.  moles/gm.  H.F.  at  C.P,  C* 


CQ* 

.007110 

94052  cal  ./mole. 

13.758  oal./molo.  °K 

%0 

.014768 

57802  - 

11.620  « 

•• 

CO 

.0104-62 

26413  - 

8.397  - 

•• 

H* 

.002478 

- 

8.030  .. 

•• 

0t 

.001754 

mm 

8.848  M 

•• 

OH 

.003107 

-  9903  -  - 

8.127  .. 

m 

H 

.001054 

-52002  .. 

4.968  .. 

— 

0 

.000791 

-59039  -  - 

4.968  .. 

-  ,  'f' 

(A 

Total 

.041530 

1666  oal./gu. 

•42714  oal^/fiole,  *K 

Nett  heat  available  for  heating  1666  -  199  ■  1467  oal«/gn« 
Heat  repaired  to  raise  temperature  to  3600*X  *  1410  oal,/gn. 


Owing  to  the  rapid  inorease  of  dissociation  with  temperature  the  effective 
specific  heat  will  be  of  the  order  of  1  oal./ga.  °K  and  so  the  surplus  of 
57  cal./ gm.  will  raise  the  gas  tenperature  try  only  50  or  60  °K.  The  xttxt  stage 

is  therefore  to  recompute  the  gas  composition  for  a  temperature  of  3660°K  using 
a  slightly  lower  value  of  y,  about  17  atmospheres,  to  allow  for  the  increased 
dissociation.  If  the  heat  balance  is  still  in  error  Ihe  correct  value  of  Tq 
and  corresponding  gas  composition  can  be  obtained  by  interpolation.  Our  next 
step  is  to  calculate  the  total  entropy,  including  the  terms  due  to  pressure  and 
to  mixing,  as  given  by  equation  (62)  and  using  table  6, 

To  calculate  conditions  at  some  stage  during  the  adiabatic  expansion  we  choose 
a  temperature  say  2000°K  and  then  calculate  the  gas  composition  and  total  pressure 
for  two  or  more  values  of  y  as  before.  For  each  value  of  y  wo  calculate  the  total 
entropy  and  then  obtain  by  interpolation  the  value  of  y,  and  the  corresponding  gas 
composition  and  total  pressure,  for  which  the  entropy  is  the  same  as  for  the 
unexpanded  gas.  I ty  repetition  at  other  temperatures,  and  interpolation  if 

necessary  we  can  then  find  the  gas  temperature  and  composition  for  any  desired 
final  pressure,  usually  one  atmosphere.  The  difference  between  the  heat  of 
reaction  and  the  heat  content  then  gives  the  surplus  enthalpy  h  converted  into 
kinetic  energy,  the  gas  velocity  being  given  by 


from  which  the  specific  impulse  v/g  may  be  calculated. 

For  the  cooler  conditions  near  the  rocket  exit  dissociation  will  be  small 
and  the  first  method  could  be  used  but  it  offers  no  great  advantage.  For 
temperatures  over  3000°K  it  is  almost  useless. 

We  could  have  shortened  the  working  slightly  by  assuming  a  constant  gas 
composition  during  the  expansion.  The  calculated  specific  impulse  would  have 
been  only  a  few  percent  lower  but  the  calculated  temperature  at  the  nozzle  exit 
would  have  been  some  500°  or  so  lower.  The  more  accurate  method  is  essential 
for  calculating  nozzle  exit  area. 
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Table  1. 


Calorim# trio  value s  and  heats  of  formation  from  the  elements 


At  300* K  and  corrected  to  lcnr  pressures. 


N.B.  The  elementary  form  of  carbon  is  taken  to  be  graphite. 


Substance 

Formula 

Mol. 

wt. 

Atomic  composition 
gn.  atcmVgq. 

Cal. rali 
oal./gm. 
oonst. 
vol. 

Heat  of  formation 
oel./gnu 

C 

H 

4 

0 

oonst. 

vol. 

oonst. 

—  Pr* 

Nitrocellulose 

ii.5£  n 

.023334 

.030681 

.008210 

.035864 

800 

668 

690 

(For  British 

12.0 

.022740 

.029333 

.008567 

.036083 

870 

636 

658 

propellants 

12.2 

.022502 

.028794 

.008709 

.036170 

898 

624 

646 

add  0,10  to  the 

12.3 

.022383 

.028525 

.008781 

.036214 

912 

618 

640 

nominal  %  N) 

12.4 

.022264 

.028255 

.008852 

.036258 

926, , 

612 

634 

12.6 

.022026 

.027716 

.008995 

.036345 

954'  * 

599 

622 

12.8 

.021789 

.027177 

.009138 

.036433 

982 

587 

609 

13.0 

.021551 

.026638 

.009280 

.036520 

1010 

574 

596 

13.1 

.021432 

.026369 

.009352 

.036564 

1024 

568 

590 

13.2 

.021313 

.026099 

.009423 

.036607 

1038 

562 

584 

13.3 

.021194 

.025832 

.009495 

.036651 

1052 

556 

578 

13.5 

.020957 

.025290 

.009637 

.036739 

1080 

543 

565 

Nitroglycerin 

CjHsNjO# 

227.09 

.013210 

.022017 

.013210 

.039631 

1750 

383 

405 

C4H,N20t 

196.12 

.020396 

.040791 

.010198 

.035692 

1050 

525 

551 

Dinitroglyool 

CaH4N,0. 

152.07 

.013152 

.026304 

.013152 

.039456 

1740 

384 

408 

Picrite 

C  H4N4O2 

104.07 

.009609 

.038434 

.038434 

.019217 

715 

189 

218 

Guanidine  nitrate 

CH4N403 

122.09 

.008191 

.049144 

.032763 

.024572 

610 

713 

745 

Oxazoide 

CtEiNtO* 

88.07 

.022710 

.045419 

.022710 

.022710 

-745 

1352 

1379 

Dinitrotolnene 

CtH*Ni04 

182.13 

.038433 

.032943 

310981 

.021962 

-135 

52 

72 

T.N.T, 

CtHsNjO* 

227.13 

.030819 

.022013 

.013208 

.026416 

500 

27 

45 

Diethyl  phthalate 

Ci  *H|  404 

222.23 

.053998 

.062997 

- 

.017999 

-1765 

782 

606 

Dibutyl  " 

Ci  aH**04 

278.34 

.057484 

.079041 

- 

.014371 

-2070 

700 

728 

Dianyl  n 

Ci  *Hi  404 

306.39 

.058749 

.084860 

•# 

.013055 

-2190 

680 

709 

Carbamite 

Ci  rUioNtO 

268.35 

.063351 

.074531 

.007453 

.003727 

-2440 

115 

141 

Methyl  central ite 

C|  jHi  «Nt  0 

240.29 

.062424 

.066585 

.008323 

-.004162 

-238O 

122 

146 

Akardite  • 

Ci,H,sNaO 

212.24 

.061251 

.056539 

.009423 

.004712 

-2140 

-34 

-13 

Mineral  jelly 

Ci  ahsa 

234.48 

.070731 

.149322 

- 

- 

-3305 

428 

473 

Special  M,J. 

- 

- 

.070791 

.148608 

• 

• 

-3260 

381 

4 26 

Diphenyls mine 

C,  «H,,N 

169.22 

.070915 

.O65OO6 

.005910 

- 

-2660 

-224 

-203 

Canphor 

C,0H,  *0 

152.23 

.065691 

.105105 

- 

.OO6569 

-2670 

442 

475 

Graphite 

C 

12.01 

.083264 

— 

- 

- 

-3365 

-51 

-21 

Water  *  ' 

HjO 

18.M6 

.111012 

_ 

.055506 

0 

3740 

3790 

Aoetone 

C9H«0 

38.08 

.051655 

.103309 

- 

.017218 

-1930 

989 

1025 

Ethyl  alcohol 

C«H*0 

46.07 

.043414 

.130242 

- 

.021707 

-1700 

1397 

1442 

Ethyl  ether 

C4H,oO 

74.12 

.053967 

.134916 

- 

.013492 

-2100 

814 

858 

T.N.T. 

CtHjNjO* 

227.13 

.O3O819 

.022013 

.013208 

.026416 

500 

27 

45 

To  try  1 

C7H5N5O. 

287.15 

.024378 

.017413 

.017413 

.027860 

930 

-46 

-28 

R.D.X. 

C3H*N404 

222.13 

.013506 

.02701? 

.027012 

.027012 

1360 

-89 

-65 

P.E.T.N, 

CsH*N40i  t 

316.15 

.015815 

.025305 

.012652 

.037957 

1535 

379 

401 

Methyl  nitrate 

CHJNDJ 

77.04 

.012980 

.038940 

.012980 

.038940 

1640 

456 

483 

Ethyl  nitrate 

CtHsNO, 

91.07 

.021962 

.054904 

.010981 

.032942 

850 

476 

505 

Nitrooe  thane 

CH3  NO* 

61.04 

.016382 

.049146 

.016382 

.032764 

1120 

421 

450 

Nitropropane 

C,H7N0* 

89.09 

.033672 

.078569 

.011224 

.022448 

-305 

448 

482 

hydrogen  peroxide 

85%,  later  15* 

- 

.066628 

- 

.058302 

748 

3180 

3207 

Nitrie  aoid 

HNOj 

63.02 

- 

.015869 

iOl 5869 

.047607 

2547 

661 

685 

Liquid  oaygen 

Os 

32.00 

- 

* 

- 

.062500 

4133 

78 

96 

Do cane 

C,  0H2X 

142.28 

.070286 

.154629 

- 

- 

-3319 

460 

506 

Diamond 

.0832ft*. 

- 

- 

- 

38 

38 

Table  2. 


Heats  of  POrmatlon  of  Explosion  Products 
In  calories  per  stole  corrected  to  loir  pressure. 


Substance 

Mol.  wt. 

-AH f 

T  -  300*K 

-  A  1 

T  ■  300*K 

3f 

T  =  0*K 

CO, 

44.010 

94052 

94052 

93969 

CO 

28.010 

26413 

26711 

27202 

H,0  liquid 

18.016 

68303 

67409 

- 

HzO  gas 

18.016 

57802 

57504 

57105 

OH 

17.008 

-9903 

-9903 

-9306 

HO 

30.008 

-21600 

-21600 

-21475 

H 

1.008 

-52002 

-51794 

-51620 

N 

14.008 

—85568 

-85270 

-85120 

0 

16.000 

-59039 

-58741 

-58590 

C04 

16.042 

17903 

17307 

15987 

HE, 

17.032 

10448 

10150 

9431 

Tabla  3. 

Rroe  energy  ftinotlcra  (HS  —  T*VT  at  cob  a  too  sphere  pywuf 

Id  oAl«/aol«  *1 


T 

00k 

H.0 

00 

& 

14 

Ot 

OH 

HO 

Sraphitt 

CBt 

Kb 

300 

400 

500 

43.605 

45.831 

47.662 

37.226 

39.514 

41.300 

40.391 

42.393 

43.947 

24.465 

26.422 

27.950 

38.859 

40.861 

42J415 

42.103 

44.111 

45.676 

36.860 

38.905 

40.484 

43.015 

45.123 

46.753 

.523 

.825 

1.146 

36.505 

38.852 

40.750 

38.185 

40.527 

42.399 

600 

700 

800 

900 

1000 

49.240 

50.638 

51.896 

53.046 

54.108 

42.774 

44.033 

45.138 

46.127 

47.025 

45.222 

46.308 

47.255 

48.098 

48.859 

29.203 

30.265 

31.186 

32.002 

32.736 

43.688 

44.768 

45.711 

46.550 

47.306 

46.967 

48.072 

49.043 

49.911 

50.697 

41.773 

42.861 

43.805 

44.638 

45.385 

48.084 

49.212 

50.197 

51.071 

51.861 

1.477 

1.810 

2.138 

2.459 

2.771 

42.389 

43.859 

45.210 

46.471 

47.662 

43.977 

45.361 

46.603 

47.736 

48.785 

1100 

1200 

1300 

1400 

1500 

55.094 

56.017 

56.883 

57.700 

58.474 

47.843 

48.613 

49.326 

49.996 

50.630 

49.555 

50.196 

50.790 

51.344 

51.864 

33.402 

34.012 

34.575 

35.100 

35.592 

47.995 

48.629 

49.218 

49.768 

50.283 

51.416 

52.079 

52.695 

53.270 

53.809 

46.063 

46.685 

47.260 

47.795 

48.296 

52.581 

53.243 

53.855 

54.426 

54.962 

3.073 

3.365 

3.647 

3.919 

4.181 

48.794 

49.873 

50.909 

51.905 

52.863 

49.764 

50.686 

51.560 

52.391 

53.184 

1600 

1700 

1800 

1900 

2000 

59.209 

59.910 

60*580 

61.221 

61.835 

51.231 

51.803 

52.350 

52.874 

53.378 

52.354 

52.818 

53.259 

53.678 

54.078 

36.055 

36.492 

36.905 

37.297 

37.671 

50.769 

51.229 

51.666 

52.082 

52.478 

54.317 

54.798 

55.254 

55.688 

56.103 

48.766 

49.210 

49.632 

50.034 

50.416 

55.466 

55.942 

56.395 

56.826 

57.237 

0 

H 

If 

2100 

2200 

2300 

2400 

2500 

62.425 

62.993 

63.540 

64.067 

64.576 

53.863 

54.331 

54.783 

55.219 

55.642 

54.460 

54.826 

55.178 

55.517 

55.843 

38.029 

38.373 

38.703 

39.021 

39.328 

52.857 

53.220 

53.569 

53.904 

54.227 

56.500 

56.880 

57.245 

57.597 

57.936 

50.781 

51.132 

51.470 

51.794 

52.106 

57.629 

58.005 

58.365 

58.711 

59.045 

43.251 

43.486 

43.711 

43.927 

44.133 

32.122 

32.353 

32.574 

32.786 

32.988 

41.3W- 

41.575 

41.796 

42.008 

42.211 

jiliit 

65.069 

65.547 

66.010 

66.459 

66.896 

56.053 

56.452 

56.839 

57.216 

57.583 

56.157 

56*461 

56.755 

57.039 

57.314 

39.625 

39.912 

40.189 

40*458 

40.719 

54.539 

54.840 

55.131 

55.413 

55.686 

58.263 

58.579 

58.885 

59.181 

59.468 

52*408 

52.701 

52.983 

53.256 

53.522 

59.366 

59.676 

59.976 

60.266 

60.548 

44.330 

44.521 

44.704 

44.882 

45.054 

33.183 

33.371 

33.551 

33.726 

33.894 

42*406 

42.594 

42.774 

42.949 

43.117 

3100 

3200 

3300 

3400 

3500 

67.320 

67.732 

68.134 

68.526 

68.907 

57.941 

58.290 

58.631 

58.964 

59.288 

57.581 

57.841 

58.094 

58.340 

58.579 

40.973 
41.220 
41.460 
41 .694 
41.922 

55.951 

56.209 

56.460 

56.704 

56.941 

59.748 
60*021 
60.28 6 
60.544 
60.794 

53.781 

54.033 

54.277 

54.514 

54.744 

60.822 
61.068 
61 .346. 
61.597 
61.841 

45.219 

45.378 

45.533 

45.683 

45.829 

34.057 

34.215 

34.368 

34.516 

34.660 

43.280 

43*438 

43.592 

43.740 

43.884 

3600 

3700 

3800 

3900 

4000 

69.279 

69.642 

69.996 

70.343 

70.683 

59.605 

59.915 

60.219 

60.517 

60.809 

58.811 

59.038 

59.260 

59.477 

59.688 

42.145 

42.363 

42.576 

42.784 

42.988 

57.172 

57.398 

57.618 

57.833 

58.043 

61.038 

61.276 

61.508 

61.735 

61.958 

54.969 

55.190 

55.406 

55.617 

55.823 

62.078 

62.309 

62.534 

62.753 

62.967 

45.971 

46.109 

46.244 

46.375 

46.504 

34.800 

34.936 

35.069 

35.198 

35.323 

44.025 

44.161 

44.294 

44.425 

44.551 

4500. 

5000 

72.273 

73.714 

62.189 

63.445 

60.674 

61.566 

43.942 

44.809 

59.024 

59.910 

63.002 

63.947 

56.790 

57.667 

63.967 

64.868 

47.098 

47.635 

35.908 

36.432 

45.142 

45.675 

Tabic  4. 

Equilibrium  Constanta 


For  partial  pressures  expressed  in  atmospheres. 


T  °X 

Ko 

log  X, 

log  Kt 

log  Xs 

log  *4 

log  K, 

log  K4 

log  XT 

log  K« 

log  Xt 

300 

400 

500 

1*118  •  10“5 
6.78  .  10-4 
.00763 

35.6293 

"ST.  8451 

£5.2130 

27.8479 

55.1207 

3^.3183 

55.6673 

55.6372 

55.1653 

£5.3458 

55.4225 

£5.2264 

24.696 

15.611 

10.060 

2.9491 

.8831 

1.5910 

"57.1881 

lt.7163 

9.2486 

600 

700 

800 

900 

1000 

.03704 

.1113 

.2486 

.4566 

.7316 

9.3560 

12.7048 

15.8685 

15.4580 

25.2164 

27.8737 

6.298 

3.5723 

1.5074 

7.8870 

5.5817 

1.7034 

2.0466 

2*5430 

2.1452 

£.8222 

5.2709 

5.4261 

2.0375 

1.2854 

0.2784 

1100 

1200 

1300 

1400 

1500 

1.065 

1.443 

1.852 

2.281 

2.718 

5.2448 

7.7032 

9.7817 

15.7740 

12.9351 

72.5465 

1.5087 

5.6ll6 

5.8518 

5.1998 

6.6340 

1.5551 

4.3312 

5.1412 

2.9774 

5.8349 

1.0870 

1.7570 

2.3200 

2.7992 

3.2118 

1600 

1700 

1800 

1900 

2000 

3.154 

3.584 

4.004 

4.411 

4.803 

5.1670 

4.5526 

£.8984 

5.2100 

£.2121 

5.7551 

7.9470 

7.3232 

5.9155 

5.1392 

7.7034 

7.3175 

5.9739 

5.6655 

2.7100 

5.6002 

5.5033 

5.4173 

5.3403 

3.5702 

3.8840 

4.1608 

4J4066 

4.6262 

2100 

2200 

2300 

2400 

2500 

5.177 

5.533 

5.870 

6.188 

6.488 

1.4924 

1.7496 

3.9848 

2.2006 

2.3993 

5*5706 

5.8968 

1.1947 

2.4678 

3.7191 

5.2453 

£.5438 

5.8584 

5.9658 

5.5469 

5.5046 

5.0512 

5.4410 

2600 

2700 

2800 

2900 

3000 

6.769 

7.033 

7.279 

7.509 

7.723 

1.5831 

1.7535 

2.9119 

1.0596 

T.1975 

2.9511 

2.1659 

2.3652 

2.5507 

2.7240 

£.8191 

2.0741 

2.3110 

2.5316 

3.7376 

5*1531 
5.4263 
it. 6799 
it.9164 
3.1372 

5.3730 

5.7503 

2.1008 

2.4275 

3.7325 

5.2382 

it.6771 

1.0020 

3.3053 

75.7469 

5.8654 

5.9570 

3100 

3200 

3300 

3400 

3500 

7.923 

8.112 

8.290 

8.455 

8.60 6 

T.3266 

T.4479 

7.5619 

1.6693 

T.7706 

2.8862 

T.0382 

T.1807 

7.3147 

1.4411 

2.9304 

2.1111 

2.2807 

2.4403 

2.5910 

2.3438 

2*5377 

2*7200 

2.8916 

2.0535 

2.0178 

2.2852 

2.5364 

2.7729 

2.9962 

2.5894 

2*8559 

2.1062 

2.3418 

2.5643 

!5.2i87 

5.7380 

3600 

3700 

3800 

3900 

4000 

8.744 

8.872 

8.994 

9.107 

9.208 

T.8663 

“.9570 

0.0430 

0.1246 

0.2019 

T.56O6 

T.6738 

T.7811 

T.8828 

T.9793 

1.7334 

2.8680 

1.9953 

1.1159 

T.2305 

I.2065 

1.3514 

2.4888 

2.6193 

1.7432 

T.2074 

T.4073 

T.5967 

1.7764 

1.9474 

1.7747 

2.9738 

1.1622 

1.3409 

1.5112 

TT.8305 

5.6459 

4500 

5000 

9.632 

9.902 

0.5389 

0.0094 

0.3864 

0.7309 

“.7259 

0.1226 

T.2820 

T.7150 

O.6865 

1.2807 

0.2465 

0.8377 

77.5394 

77.3108 

5.5780 

5.5264 

N.B.  For  explanation  see  next  page 


Explanation  of  Table 


*•  -Mf?m  -  rntm 

k,  =  p(h)/v1p(iOJ 

K*  =  p(OH)  /{pC^JJ/pCBeO) 

K3  =  p(ND)  .  p(H«)/p(H«0)  ^p(N2)| 

K4  =  p(n)/v1p(n2)| 

K5  =  p(0)  .  pCHjVpCHaO) 

K«  =  p(02)  .  jpC^VpCHaO)!2 
Kt  =  PCCHa)  •  p(HaO)/p(CO)  {p(Ha)]s 
K.  =  p(NIfc)/Vh>(N2)]  •  {P(H*)]S/2 

K9  =  jp(CO) j2/p(C02 )  in  presence  of  graphite. 


The  gases  are  assumed  perfect  in  all  oases. 


Table  5. 

_ 

Mean  Specific  Heats  Oy  between  300  and  T  °K  at  Constant  Voluiae 

In  oal./mole  °K 


T  #K 

GOt 

H*0 

00 

H> 

Nt 

o> 

OH 

NO 

Graphite 

ch4 

NHj 

0 

300 

400 

500 

5.527 

6.904 

7.412 

7.853 

5.950 

6.042 

6.115 

6.212 

4,964. 

4.978 

4.996 

5.035 

4.801 

4.908 

4.958 

4.978 

4.963 

4.974 

4.986 

5.0J3 

4.955 

5.032 

5.110 

5.220 

5.077 

5.153 

5.112 

5.091 

5.377 

5.158 

5.159 
5.195 

.851 

2.083 

2.473 

2.828 

6.052 

6.586 

7.134 

7.799 

6.080 

6.521 

6.843 

7.208 

600 

700 

800 

900 

1000 

8.238 

8.576 

8.876 

9.145 

9.385 

6.328 

6.452 

6.586 

6.725 

6.868 

5.093 

5.164 

5.241 

5.321 

5.400 

4.987 

5,000 

5.015 

5.033 

5.055 

5.056 

5.113 

5.180 

5.251 

5.323 

5.334 

5.448 

5.556 

5.656 

5.747 

5.081 

5.080 

5.090 

5.108 

5.131 

5.265 

5.346 

5.429 

5.513 

5.594 

3.143 

3.417 

3.653 

3.855 

4.027 

8.497 

9.192 

9.866 

10.510 

11.122 

7.581 

7.951 

8.315 

8.671 

9.016 

1100 

1200 

1300 

1400 

1500 

9.601 

9.796 

9.972 

10.132 

10.278 

7.013 

7.159 

7.304 

7.448 

7.589 

5.477 

5.549 

5.617 

5.681 

5.741 

5.082 

5.114 

5.150 

5.188 

5.228 

5.395 

5.465 

5.532 

5.595 

5.654 

5.829 

5.903 

5.971 

6.033 

6.091 

5.160 

5.196 

5.240 

5.288 

5.336 

5.670 

5.741 

5.807 

5.868 

5.925 

4.175 

4.306 

4.425 

4.533 

4.632 

11.700 

12.243 

12.752 

13.228 

13.672 

9.349 

9.669 

9.975 

10.268 

10.546 

1600 

1700 

1800 

1900 

2000 

10.4.11 

10.533 

10.645 

10.748 

10.844 

7.726 

7.860 

7.989 

8.115 

8.236 

5.797 

5.849 

5.897 

5.942 

5.984 

5.270 

5.313 

5.357 

5.401 

5.445 

5.710 

5.763 

5.813 

5.860 

5.903 

6.144 

6.194 

6.241 

6.286 

6.328 

5.383 

5.430 

5.475 

5.519 

5.562 

5.977 

6.025 

6.068 

6.108 

6.146 

14.089 

14.477 

14.840 

15.180 

15.497 

10.811 

11.062 

11.299 

11.524 

11.737 

2100 

2200 

2300 

2400 

2500 

10.933 

11.015 

11.092 

11.164 

11.231 

8.352 

8.464 

8.572 

8.674 

8.773 

6.024 

6.061 

6.096 

6.128 

6.159 

5.489 

5.532 

5.574 

5.615 

5.656 

5.943 

5.981 

6.017 

6.051 

6.O83 

6.369 

6.409 

6.447 

6.484 

6.520 

5.604 

5.645 

5.686 

5.727 

5.767 

6.182 

6.216 

6.248 

6.278 

6.306 

16.819 

1 2.644 

2600 

2700 

2800 

2900 

3000 

11.295 

11.355 

11.411 

11.464 

11.514 

8.868 

8.959 

9.046 

9.130 

9.210 

6.189 

6.217 

6.243 

6.267 

6.290 

5.696 

5.735 

5.773 

5.810 

5.846 

6.114 

6.143 

6.171 

6.197 

6.221 

6.555 

6.589 

6.622 

6.654 

6.685 

5.809 

5.846 

5.884 

5.920 

5.955 

6.332 

6.357 

6.381 

6.404 

6.425 

6.0 

17.807 

13.343 

3100 

3200 

3300 

3400 

3500 

11.562 
11 .608 
11.651 
11.693 
11.733 

9.288 

9.362 

9.434 

9.503 

9.569 

6.312 

6.334 

6.355 

6.374 

6.392 

5.881 

5.915 

5.948 

5.981 

6.013 

6.245 

6.268 

6.290 

6.310 

6.329 

6.716 

6.746 

6.776 

6.805 

6.833 

5.988 

6.021 

6.053 

6.O83 

6.112 

6.445 

6.465 

6.484 

6.503 

6.521 

18.570 

13.893 

3600 

3700 

3800 

3900 

4000 

11.771 

11.807 

11.842 

11.876 

11.908 

9.633 

9.694 

9.753 

9.810 

9.864 

6.410 

6.427 

6.444 

6.460 

6.475 

6.043 

6.072 

6.101 

6.129 

6.156 

6.348 

6.366 

6.383 

6.400 

6.416 

6.861 

6.889 

6.916 

6.942 

6.968 

6.140 

6.168 

6.195 

6.221 

6.246 

6.536 

6.552 

6.567 

6.582 

6.597 

19.175 

14.337 

4500 

5000 

12.053 

12.177 

10.110 

10.310 

6.542 

6.602 

6.285 

6.398 

6.490 

6.548 

7.084 

7.191 

6.372 

6.482 

6.668 

6.730 

19.667 

20.075 

14.701 

15.007 

For  monatomic  gases  take  Cy  =  2.981. 
To  obtain  Cp  add  1 .987. 


Entropies  at  One  Atmosphere  Pressure 

In  oal./mole  °K 


T  °K 

COa 

HtO 

CO 

H. 

N* 

Oa 

OH 

NO 

Graphite 

CH4 

NHj 

300 

400 

500 

51.119 

53.816 

56.107 

45.162 

47.492 

49.342 

47.342 

49.352 

50.927 

31.253 

33.250 

34.809 

45.809 

47.818 

49.385 

49.045 

51.092 

52.724 

43.924 

45.978 

47.554 

50.379 

52.433 

54.044 

1.374 

2.081 

2.788 

45.078 

48.600 

51.917 

46.252 

48.785 

50.917 

600 

700 

800 

900 

1000 

58.110 

59.893 

61.503 

62.972 

64.323 

50.900 

52.257 

53.472 

54.581 

55.605 

52.238 

53.373 

54.379 

55.287 

56.115 

36.084 

37.167 

38.108 

38.945 

39.702 

50.685 

51.805 

52.797 

53.692 

54.509 

54.099 

55.296 

56.361 

57.321 

58.194 

48.839 

49.927 

50.878 

51.723 

52.487 

55.392 

56.558 

57.593 

58.526 

59.377 

3.474 

4.127 

4.740 

5.313 

57!732 

6O.364 

62.743 

64.938 

52.795 

54.497 

56.067 

57.531 

58.907 

1100 

1200 

1300 

1400 

1500 

65.572 

66.733 

67.817 

68.833 

69.789 

56.559 

57.456 

58.304 

59.110 

59.878 

56.878 

57.586 

58.245 

58.860 

59.437 

40.395 
41 .035 
41 .631 
42.192 
42.722 

55.260 

55.956 

56.606 

57.215 

57.786 

58.994 
59.732 
60.418 
61 .059 
61 .660 

53.188 

53.838 

54.450 

55.025 

55.567 

60.158 

60.880 

61.550 

62.176 

62.764 

6.342 

6.8O7 

7.247 

7.663 

8.057 

66.963 

68.835 

70.570 

72.182 

73.682 

6^*08  v 

61.445 

62.624 

63.749 

1600 

1700 

1800 

1900 

2000 

70.692 

71.547 

72.359 

73.132 

73.869 

60.611 

61.313 

61.986 

62.634 

63.258 

59.982 
60.498 
60.988 
61 .453 
61.896 

43.224 

43.702 

44.157 

44.591 

45.007 

58.326 

58.838 

59.325 

59.788 

60.228 

62.226 

62.761 

63.268 

63.751 

64.212 

56.079 

56.565 

57.028 

57.470 

57.892 

63.318 

63.840 

64.335 

64.806 

65.255 

0 

H 

N 

2100 

2200 

2300 

2400 

2500 

74.573 

75.247 

75.893 

76.513 

77.110 

63.859 

64.439 

64.999 

65.540 

66.064 

62.320 

62.726 

63.114 

63.486 

63.845 

45.407 

45.792 

46.163 

46.521 

4-6.868 

60.647 
61.049 
61 .435 
61. 806 
62.163 

64.654 

65.078 

65.485 

65.877 

66.255 

58.297 

58.687 

59.064 

59.427 

59.778 

65.684 

66.094 

66.487 

66.864 

67.226 

48.219 

48.454 

48.679 

48.895 

49.101 

37.090 

37.321 

37.542 

37.754 

37.956 

46.312 

46.543 

46.764 

46.976 

47.179 

2600 

2700 

2800 

2900 

3000 

77.686 

78.242 

78.779 

79.298 

79.800 

66.571 

67.063 

67.540 

68.004 

68.455 

64.191 

64.525 

64.047 

65.158 

65.459 

47.204 

47.530 

47.845 

48.151 

48.448 

62.507 

62.839 

63.16O 

63.470 

63.770 

66.621 

66.974 

67.315 

67.646 

67.968 

60.119 
60.448 
60.767 
61.076 
61 .376 

67.575 

67.911 

68.236 

68.551 

68.856 

49.298 

49.489 

49.672 

49.850 

50.022 

38.151 

38.339 

38.518 

38.694 

38.862 

47.374 
47.562 
47.742 
47.91 7 
48.085 

3100 

3200 

3300 

3400 

3500 

80.286 

80.758 

81.216 

81.662 

82.096 

68.894 

69.321 

69.736 

70.140 

70.533 

65.751 

66.034 

66.309 

66.576 

66.835 

48.737 

49.018 

49.291 

49.558 

49.818 

64.O6O 

64.342 

64.616 

64.882 

65.141 

68.282 

68.587 

68.884 

69.172 

69.453 

61 .668 
61.952 
62.228 
62.495 
62.755 

69.151 

69.438 

69.717 

69.988 

70.251 

50.187 

50.346 

50.501 

50.651 

50.797 

39.025 

39.183 

39.336 

39.484 

39.628 

48.248 

48.406 

48.560 

48.708 

48.852 

3600 

3700 

3800 

3900 

4000 

82.517 

82.927 

83.326 

83.716 

84.099 

70.917 

71.292 

71.659 

72.018 

72.368 

67.087 

67.333 

67.573 

67.8O8 

68.037 

50.072 
50.319 
50.561 
50.798 
51 .030 

65.392 

65.637 

65.876 

66.110 

66.338 

69.727 

69.995 

70.256 

70.511 

70.762 

63.009 

63.257 

63.500 

63.737 

63.969 

70.506 

70.753 

70.994 

71.230 

71.460 

50.939 

51.077 

51.212 

51.343 

51.472 

39.768 

39.904 

40.037 

40.166 

40.291 

48.993 

49.129 

49.262 

49.393 

49.519 

4500 

5000 

85.878 

87.479 

74.008 

75.481 

69.100 

70.056 

52.115 

53.099 

67.402 

68.351 

71 .931 
72.991 

65.O63 

66.052 

72.536 

73.504 

52.066 

52.603 

40.876 

41.400 

50.110 

50.643 

Table  7 


Specific  Heats  Cy  at  Constant  Volume 
In  cal ./mole  °K 


T  °K 

COa 

HaO 

CO 

H* 

Na 

Oa 

OH 

NO 

Graphite 

CH4 

NHs 

300 

400 

500 

6.904 

7.885 

8.672 

6.042 

6.202 

6.430 

4.978 

5.026 

5.133 

4. 908 
4.987 
5.006 

4.974 

5.004 

5.083 

5.032 

5.207 

5.432 

5.153 

5.086 

5.061 

5.158 

5.180 

5.305 

2.083 

2.851 

3.506 

6.586 

7.761 

9.179 

6.521 

7.196 

7.949 

600 

700 

800 

900 

1000 

9.316 

9.850 

10.296 

10.668 

10.979 

6.695 

6.974 

7.267 

7.569 

7.873 

5.289 

5.464 

5.637 

5.800 

5.945 

5.021 

5.048 

5.091 

5.152 

5.230 

5.210 

5.364 

5.526 

5.684 

5.829 

5.683 

5.898 

6.077 

6.225 

6.348 

5.064 

5.099 

5.162 

5.247 

5.346 

5.486 

5.673 

5.848 

6.006 

6.140 

4.030 

4.440 

4.750 

4.980 

5.140 

10.599 

11.936 

13.166 

14.284 

15.288 

8.700 

9.423 

IO.115 

10.772 

11.390 

1100 

1200 

1300 

1400 

1500 

11.241 
11 .462 
11.649 
11.808 
11.945 

8.175 

8.468 

8.749 

9.015 

9.263 

6.071 

6.181 

6.277 

6.360 

6.432 

5.320 

5.417 

5.518 

5.622 

5.726 

5.960 

6.076 

6.177 

6.266 

6.344 

6.453 

6.543 

6.621 

6.690 

6.752 

5.454 

5.566 

5.676 

5.784 

5.888 

6.255 

6.354 

6.438 

6.509 

6.571 

5.280 

5.420 

5.560 

5.680 

5.760 

16.182 

16.974 

17.673 

18.290 

18.832 

11.965 

12.493 

12.974 

13.410 

13.804 

1600 

1700 

1800 

1900 

2000 

12.064 

12.167 

12.257 

12.337 

12.408 

9.494 

9.707 

9.903 

10.084 

10.250 

6.494 

6.548 

6.597 

6.640 

6.678 

5.827 

5.924 

6.017 

6.IO5 

6.188 

6.413 

6.473 

6.526 

6.573 

6.615 

6.811 

6.868 

6.925 

6.981 

7.037 

5.986 

6.078 

6.163 

6.243 

6.318 

6.625 

6.673 

6.715 

6.753 

6.787 

2100 
2200 
2300 
2400 
•  2500 

12.471 

12.528 

12.580 

12.627 

12.670 

10.403 

10.544 

10.675 

10.795 

IO.905 

6.712 

6.743 

6.771 

6.797 

6.821 

6.267 

6.341 

6.411 

6.477 

6.539 

6.653 

6.687 

6.718 

6.746 

6.772 

7.093 

7.150 

7.208 

7.265 

7.318 

6.388 

6.453 

6.514 

6.571 

6.625 

6.816 

6.842 

6.866 

6.888 

6.909 

Fo 

e 

C 

To 

a 

r  monato 
ases  tak 
v  =  2*98 
obtain 

aa  1 .987 

mio 

e 

1. 

Op 

2600 

2700 

2800 

2900 

3000 

12.710 

12.747 

12.781 

12.813 

12.843 

11.007 
11 .102 
11.189 
11.269 
11.344 

6.843 

6.863 

6.881 

6.897 

6.912 

6.598 

6.654 

6.707 

6.757 

6.804 

6.796 

6.818 

6.839 

6.858 

6.875 

7.365 

7.408 

7.450 

7.491 

7.531 

6.675 

6.723 

6.768 

6.811 

6.852 

6.928 

6.946 

6.962 

6.978 

6.993 

3100 

3200 

3300 

3400 

3500 

12.871 

12.898 

12.924 

12.949 

12.972 

11.413 

11.478 

11.538 

11.593 

11.645 

6.926 

6.940 

6.953 

6.965 

6.976 

6.848 

6.890 

6.930 

6.969 

7.006 

6.891 

6.906 

6.920 

6.934 

6.947 

7.571 

7.611 

7.651 

7.687 

7.724 

6.891 

6.928 

6.963 

6.997 

7.030 

7.008 

7.022 

7.035 

7.048 

7.061 

3600 

3700 

3800 

3900 

4000 

12.995 

13.017 

13.038 

13.059 

13.079 

11.693 

11.737 

11.777 

11.814 

11.847 

6.987 

6.998 

7.008 

7.018 

7.028 

7.041 

7.074 

7.105 

7.135 

7.164 

6.959 

6.970 

6.981 

6.992 

7.002 

7.760 

7.795 

7.829 

7.861 

7.892 

7.06i 

7.091 

7.120 

7.148 

7.174 

7.074 

7.086 

7.097 

7.109 

7.120 

4500 

5000 

13.174 

13.264 

11.974 

11.986 

7.072 

7.112 

7.295 

7.402 

7.049 

7.089 

8.016 

8.128 

7.343 

7.518 

7.172 

7.221 
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